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Abstract 

We present a comprehensive analysis of the prescription we recently put forward 
for the computation of real-time correlation functions using gauge/gravity duality. The 
prescription is valid for any holographic supergravity background and it naturally maps 
initial and final data in the bulk to initial and final states or density matrices in the 
field theory. We show in detail how the technique of holographic renormalization can 
be applied in this setting and wc provide numerous illustrative examples, including the 
computation of time-ordered, Wightman and retarded 2-point functions in Poincarc and 
global coordinates, thermal correlators and higher-point functions. 
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1 Introduction 



Since the advent of the AdS/CFT correspondence [UEIIH] considerable amount of work has 
been devoted to developing holographic dualities leading to a very precise understanding 
of the holographic dictionary in Euclidean signature, see [U El [6] for reviews. In this 
paper, continuing our recent work [7], we aim at developing a real-time prescription that is 
applicable at the same level of generality as the corresponding Euclidean prescription. More 
precisely, we would like to have a setup that is valid for all QFTs that have a holographic dual 
and is applicable for the holographic computation of n-point functions of gauge invariant 
operators in non-trivial states. Furthermore, this prescription should be fully holographic, 
i.e. only boundary data and regularity in the interior should be needed for the computation, 
and within the supergravity approximation all information should be encoded in classical 
bulk dynamics. Such a prescription is an integral part of the definition of the holographic 
correspondence and as such it is important on general grounds. Furthermore, there is a 
wide range of applications for such a general real-time prescription. To mention a few: 
one would like to understand better holography for time dependent backgrounds, to have 
a holographic description of non-equilibrium QFT and to be able to compute correlators 
in non-trivial states. Such a development is also becoming urgent as potential current and 
future applications of holography to modelling of the quark-gluon plasma in RHIC and LHC 
require real-time techniques, see [8] for a review. 

From a more theoretical perspective, one would like to understand better the interplay 
between causality and holography. Since bulk and boundary lightcones are different, it is 
not a priori clear that a bulk computation will produce the correct causal structure for 
boundary correlators, for example the correct ie insertions. Conversely, one can ask how 
the bulk causal structure emerges from boundary correlators. A related question is to 
understand how black hole horizons are encoded in boundary correlators. More generally 
one would like to study holographically the process of gravitational collapse. All of these 
applications require a formalism that can handle the general case, rather than being tied 
up to particularities of specific examples. 

In many applications, the Euclidean holographic techniques for obtaining correlators are 
sufficient, since one can often analytically continue from Lorentzian to Euclidean signature. 
While such a Wick rotation is often the most direct way of arriving at the result, there 
are also many cases where the analytic continuation is technically difficult, even though 
possible in principle. For example, in the case of a thermal correlator one would need to 
continue from a discrete set of Matsubara frequencies which is technically not very easy. 
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More importantly, the Wick rotation obscures the bulk and boundary dynamics and none 
of the questions raised above can be answered in this setting. 

There have been several earlier works discussing holography in Lorentzian signature, 
including P[ini[IIl[l2l[l3l[ia[l5l[l6l[I71[l8l[l9]. One set of these papers is based on 
semi-classical bulk quantization around a classical bulk solution. For example, a case often 
discussed is that of the quantization of a free scalar field in AdS and the computation of 
the associated boundary 2-point function. Such results are clearly difficult to extend to 
cases where bulk interactions are essential because of the difficulty in quantizing the bulk 
gravitational theory. For example, higher point functions, correlators of the stress energy 
tensor and holographic RG flows are outside the remit of these works. Moreover, for the 
computation of the correlators in the large N and strong 't Hooft limit, one should not 
have to consider the quantization of the bulk theory at all ~ classical bulk dynamics should 
suffice. 

A Lorentzian prescription that has been used widely in the literature is that of Son and 
Starinets [12] . This prescription leads to the computation of retarded 2-point functions and 
is based on imposing specific boundary conditions in the interior of the spacetime. More 
specifically, this prescription assumes that one deals with a spacetime with a horizon and 
imposes incoming-wave boundary conditions at the horizon. It was later shown in [13] 
that these conditions are related to boundary conditions discussed earlier in the black hole 
literature [20^ [2T] and that (with such boundary conditions understood) the prescription 
follows from taking functional derivatives of the on-shell action, although the authors did 
not take into account contributions to the on-shell action from timelike infinity, which are 
generically non-zero. 

This prescription leads to correct results (provided the infinities have been subtracted 
correctly, see below). It is somewhat unsatisfactory, however, from the holographic point 
of view, at least in view of the general applications we have in mind. For holography one 
would want to have all information encoded in boundary data, so that boundary data is 
sufficient on its own to reconstruct the bulk dynamics. The prescription in [121 113j . on the 
other hand, presumes the existence of a horizon and uses specific behavior of the bulk fields 
there. Furthermore, as mentioned above there are additional contributions in the on-shell 
action from initial and final surfaces within the setup of [12^ [T3] . We will see that these 
additional terms cancel in our prescription and we will also derive the behavior of the fields 
at the horizon used in [12^ [T3] from a fully holographic prescription. 

Another issue that has not been discussed adequately in the past is that of renormal- 
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ization. First, there are the infinities due to the non-compactness of the bulk spacetime 
in the radial direction. These are the same infinities appearing in the Euclidean setup 
and can be dealt with by a straightforward adaptation of the Euclidean discussion, namely 
by introducing boundary counterterms etc., see [6] for a review. Such counterterms not 
only remove infinities but also in general affect the finite part of correlators, even in such 
simple cases as the 2-point functions of scalar operators. For example, the naive computa- 
tion of 2-point functions of scalar operators of scaling dimension A 7^ d, where the infinite 
terms are simply dropped, leads to results that are inconsistent with Ward identities, see 
[22I [23l [21} [6l [251 [26] for examples and discussions of this point. A second issue that is 
specific to the Lorentzian setup is that there may also be new infinities due to the non- 
compactness in the time direction and one would also have to understand how to deal with 
those. 

The main difference between the Euclidean and the Lorentzian cases is that in the latter 
case one also has to specify initial and final conditions for the bulk fields. It has long been 
appreciated [9] that these conditions should be related to a choice of in- and out-state in 
the Lorentzian boundary QFT, but the precise relation was given only recently in [7], which 
we now review. The starting point in [7] was the fact that in QFT the initial and final 
conditions can be implemented using a time contour in a complex time plane. For example, 
one can compute the expectation values in a non-trivial state in the 'in-in' formalism by 
choosing a closed time contour that starts from the operator that creates the state, runs 
along the real time axis and returns to the operator [271 [28l [29| I30j . Thermal correlators 
can be obtained by having the contour run also along the imaginary axis. 

Since the gauge/gravity duality is believed to be an exact equivalence, one should be 
able to holographically compute correlators in non-trivial states and the choice of a contour 
in the complex time plane should be refiected on the gravitational side too. The prescription 
of [7] is therefore to start from the QFT contour and 'fill it in' with a bulk manifold. Real 
segments of the contour are associated with Lorentzian spacetimes, and purely imaginary 
segments with Euclidean solutionis]. The Euclidean bulk solution which is associated with 
the initial state on the QFT side can also be thought of as providing a Hartle-Hawking wave 
function [3T] for the bulk theory. Thus our prescription is not only QFT inspired but also in 
line with standard considerations on wave functions in quantum gravity, see also [11^ [T6] for 
related discussions. There has been considerable discussion in the literature over the choice 
^With Euclidean solutions we mean a solution of the field equations after Wick rotation to positive 
definite signature. If the solution is real, it should be more properly called 'Riemannian'. We will see in 
examples, however, that the 'Euclidean' solutions can also be complex. 
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of contour in the Euclidean path integrals and the reality conditions of the the semi-classical 
saddle point evaluation, see for example [32]. In our case, the bulk reality conditions are 
dictated by the boundary theory and, in particular, for a generic complex boundary contour 
the bulk manifold would have a complex metric (but in all cases the boundary correlators 
would satisfy standard reality conditions). 

In this paper we extend and further develop the framework discussed in [7] by presenting 
a comprehensive discussion of all issues involved. Although the discussion below is meant 
to be self-contained, to fix ideas it may be helpful to read the concrete example presented in 
[7] first. The organization of this paper is the following. In the next section we present the 
general prescription in detail. Holographic renormalization is discussed in section [3l This 
section contains important but rather technical results and can be skipped on a first reading. 
In section HI we discuss a range of different examples that complement the example of [7] ■ In 
particular, we discuss the holographic computation of a Wightman function, compute a real- 
time two-point function in thermal AdS, in eternal BTZ and in AdS in Poincare coordinates. 
We also illustrate how to compute higher-point functions and discuss the prescription for 
rotating black holes. Finally, we summarize some relevant background QFT material in the 
appendix. 

2 Real-time prescription 

In this section we discuss the real-time gauge/gravity prescription in detail. We begin with a 
discussion of basic QFT results that will be the springboard for the holographic prescription 
in subsection 12.21 

2.1 QFT preliminaries 

Consider a field configuration with initial condition (j)-~{x) at t = —T and final condition 
4>+{x) at t = T. The path integral with fields constrained to satisfy these conditions pro- 
duces the transition amplitude {4>+, T\(j)^, —T). If we are interested in vacuum amplitudes 
we should multiply this expression by the vacuum wave functions (0|(/)+,T) and ((/>_, —T\0) 
and integrate over (p^ , . The insertion of these wave functions is equivalent to extending 
the fields in the path integral to live along a contour in the complex time plane as sketched 
in Fig. [1^. Indeed, the infinite vertical segment starting at —T corresponds to a transition 
amplitude lim^^oo (</>_, —T|e~^^|^') for some state |^'), which is however irrelevant since 
taking the limit projects it onto the vacuum wave function (<^_, — T|0). Similarly, we obtain 
(0|</)+,T) from the vertical segment starting at t = T. As is reviewed in appendix |Al these 
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wave function insertions ultimately lead to the ie factors in the Feynman propagators. 

In this discussion, we used the Euclidean path integral to create the vacuum state which 
is then fed into the Lorentzian path integral as the initial and final state. More generally, one 
can use the Euclidean path integral to generate other states that can serve as initial/final 
states for the Lorentzian path integral. In the context of a conformal field theory on 
the relation between Euclidean path integrals and states is the basis for the operator-state 
correspondence: inserting a local operator O, say at the origin of MJ^, and then performing 
the path integral over the interior of the sphere S'^~^ that surrounds the origin results in the 
corresponding quantum state l^'o) on S'^~^. In particular, the vacuum state is generated 
by inserting the identity operator. 

Suppose now that we are interested in computing real-time correlation functions in a 
given initial state, {^\Oi{xi) ■ ■ ■ One can do this by using a closed time contour 

[271 l28l [29l [30] . as sketched in Fig. [Da. In the figure, the vertical pieces C^Ca represent 
Euclidean path integrals, with the crosses representing operator insertions. As described 
in the previous paragraph, these segments create the chosen initial state |^). We then 
evolve this state forward and backward in time following the horizontal segments Ci and 
C2. To compute real-time expectation values, one may insert operators either in Ci or C2. 
For example, one can have 2-point functions with operators inserted both in Ci, or one in 
Ci and one in C2, or both in C2. This leads to a 2 x 2 matrix of 2-point functions, the 
Schwinger-Keldysh propagators, discussed in more detail in appendix lAl 

For real-time thermal correlators one can use the closed time path contour in Fig.[T]:. The 
vertical segment now represents the thermal density matrix, p = exp(—PH), with /? = 1/T 
and H the Hamiltonian. The circles indicate points that should be identified and reflect 
the fact that thermal correlators satisfy appropriate periodicity conditions in imaginary 
time (bosonic/fermionic correlators are periodic/antiperiodic). As in the discussion in the 
previous paragraph, one can insert operators at any point in the horizontal segments. Other 
density matrices, for example a thermal density matrix with chemical potentials, may be 
obtained in a similar manner. 

For all of the contours above, one can write a generating functional of correlation func- 
tions of gauge invariant operators in non-trivial states with the following path integral 
representation: 



1) 



Here if denotes collectively all QFT fields, cpj^^^ are sources that couple to gauge invariant 
operators , 5(o)jj is the spacetime metric (and also the source for the stress energy tensor 
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Figure 1: (a) Vacuum-to-vacuum contour, (b) In-in contour, (c) Real-time thermal contour. 
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Tij). The path integral is performed for fields living on the contour C in the complex 
time plane. Therefore, we think of t as a complex coordinate and dt is then a contour 
integral. By the usual 'slicing' arguments in deriving a path integral, one obtains contour- 
time-ordered correlators after functionally differentiating w.r.t. sources. 

Let us exemplify ()2.1.ip using the contour in Fig. [it. In this case the contour has four 
segments and the path integral splits into four corresponding parts. The segments CqjC^ 
are associated with Euclidean path integrals, as discussed above. In these segments, we can 
parametrize the contour using t = —it with r a real coordinate along C. This substitution 
leads to the usual signs in the Euclidean path integral. The source terms in these segments 
are related to the choice of initial and final states. The segments Ci and C2 form a closed 
time path. Let us parametrize this path using a contour time coordinate tc that increases 
monotonically along the contour. In the segment Ci we can simply set t = tc, where now tc 
ranges from to T (where T may be 00), and we can integrate along C2 using t = 2T — tc, 
with T < tc < 2T. (Notice that dt = —dtc on C2, giving rise to an important extra sign 
for the action on this segment.) The source terms in these segments are the usual sources, 
which upon functional differentiation lead to contour-time-ordered n-point functions. For 
operator insertions on Ci and C2 contour-time-ordering coincides with tc-time-ordering, see 
appendix |A] for further discussion. 

Notice that in quantum field theory, one may consider deforming the contour in the 
complex time plane in any direction. In fact, one may deform it into any other direction in 
complex coordinate space. Such deformations are allowed, provided the contour does not 
run upward in the complex time plane (so that the path integral converges) and similar 
restrictions apply for deformations in complex coordinate space. In general, the 'metric' 
along such a deformed contour would be complex, which may lead to a complex bulk metric 
as well. We will consider such an example when we analyze the rotating BTZ black hole in 
section 14.51 
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2.2 Prescription 

We are now ready to present the real-time gauge/gravity prescription. Our starting point 
is that the contour dependence we discussed in the previous section should be reflected 
in the bulk string theory, and in the low energy approximation it should be part of the 
supergravity description. Within the saddle-point approximation, our prescription is to 
associate supergravity solutions with QFT contours, or, more figuratively, to 'fill in' the 
QFT contour with a bulk solution. We have sketched several examples of such a construction 
in Fig. m on page [36l 

One can think of the field theory contour C as a d-dimensional subspace of a com- 
plexified boundary spacetime. In most cases, as we saw above, this would be a line in the 
complexified time plane times a real space, M x S'^"^. The bulk solution should have C as its 
conformal boundary and the bulk fields <I>^ should satisfy boundary conditions parametrized 
by fields 0|q^(x) living on C. This means that horizontal segments of C will be filled in 
with Lorentzian solutions, while vertical segments will be filled in with Euclidean solutions. 
These segments are then glued together along bulk hypersurfaces that end on the corners 
of the contour. The total manifold consisting of all these segments is denoted by Mq and it 
has a metric whose signature jumps at the 'corner' hypersurfaces where a vertical segment 
meets with a horizontal one. Below we show how appropriate matching conditions control 
the behavior of the fields at these hypersurfaces. 

Note that the bulk manifold is not necessarily of the form R x X'^ with dX'^ = S'^"^. 
Instead, we can have more general bulk solutions that may 'interpolate' between various 
parts of the contour. An important example is the eternal BTZ black hole we consider 
below. 

Given such a solution Mc that fills in the entire field theory contour C, the next step in 
the prescription is to compute the corresponding on-shell supergravity action. This action 
is then identified with the generating functional of correlators in non-trivial states discussed 
in the previous subsection, 

ZQFT[<t^U-M = exp j^^ d''^'xV^CZi:''^'\ct>[,^]^ . (2.2.1) 

Vertical segments of the contour involve the Euclidean on-shell action and horizontal seg- 
ments the Lorentzian on-shell action, with factors of i and signs becoming standard when 
one passes from the complex coordinate t to the corresponding real contour time variable 
r or tc- Notice that this discussion does not require that tc and r extend globally on M, 
as the asymptotic analysis suffices to fix all signs. The sources (plr^, that are localized in 
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the conformal boundary of the Euchdean part of the solution are associated with the initial 
and final state, whereas sources on the conformal boundary of the Lorentzian solution lead 
to n-point functions upon differentiation. Note that (|2.2.1|) is a bare relation, as both sides 
diverge. The holographic renormalization needed to render the on-shell supergravity action 
finite will be described in section [3l 

2.2.1 Corners 

Piecewise straight contours have corners, where either a horizontal and a vertical segment 
meet or two horizontal segments join. These corners extend to hypersurfaces S in the bulk. 
The signature of the metric changes at the hypersurface corresponding to a corner of a 
horizontal and a vertical segment, but otherwise it remains unchanged. Modulo subtleties 
at the boundary of S, which we discuss in the next section, we impose the following two 
matching conditions at S: 

1. We impose continuity of the fields across S. That is, we require the induced metric, 
the values of the scalars, and induced values of the other fields to be continuous; 

2. If the contour passes from a segment M_ to Af_|_, then we impose appropriate conti- 
nuity of the conjugate momenta across S: 

7r_ = r/7r+ , (2.2.2) 

where Tr± denote collectively the conjugate momenta of all fields on the two sides M± of S 
(defined using tc or r as the time coordinate), and rj = —i when we consider a Euclidean to 
a Lorentzian corner like for example from Cq to Ci in Fig. [T]3, whereas rj = —1 if we have a 
(non-trivial) Lorentzian to Lorentzian corner as from Ci to C2 in Fig. [T]3. In all cases, the 
matching condition is equivalent to 

7r+ = 7r_ , (2.2.3) 

where tc is defined using the complex time variable t. In other words, if we use analytic 
continuation of the fields in the complex t coordinate to smooth out the corner by bending 
the contour, then the matching conditions dictate that the solution would be at least C^. 
In section [J] we illustrate with examples how these matching conditions determine the bulk 
solution for a given contour. 

The on-shell supergravity action can be regarded as the saddle point approximation of 
the 'bulk path integral' and the matching conditions can also be justified in the same way. 
Recall that a path integral for fields living on a certain manifold can always be split in two 
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by cutting the manifold in two halves and imposing boundary conditions for the fields on 
the cut surface. Afterwards, one can glue the pieces back together by imposing the same 
boundary condition on either side and then integrate over these boundary conditions. 

The saddle-point approximation can similarly be performed in steps. After cutting the 
manifold, one first finds a saddle-point approximation on either side with arbitrary initial 
data at the cut surface. This replaces the partial path integrals on either side by an on-shell 
action which in particular depends on the initial data. Then, one imposes continuity of 
the initial data, which is the first matching condition, and performs a second saddle-point 
approximation with respect to the initial data. Since the first variation of an on-shell action 
with respect to boundary data yields the conjugate momentum, this second saddle-point 
precisely yields (|2.2.2p . The matching conditions should then be viewed as an equation 
determining the initial data. One may verify that the signs come out right, too. This 
formalism presupposes a two-derivative action, but higher-derivative terms can be dealt 
with perturbatively. 

3 Holographic renormalization 

The fundamental holographic relation (j2.2.ip is a bare relation because both sides are di- 
vergent: there are UV divergences on the QFT side and IR infinite volume divergences 
on the gravitational side. So appropriate renormalization is needed to make this relation 
well-defined. In this section, we will show that the procedure of holographic renormaliza- 
tion for the spaces under consideration is a priori more complicated, but that none of these 
complications enter in the final result. Therefore, the formulas presented in for example [6] 
for the holographically computed correlation functions remain valid in the context of our 
real-time prescription as well. As the precise derivation of this result is not essential for the 
rest of the paper, the reader may wish to skip this section on a first reading and proceed 
directly to the examples of section HI 

The holographic renormalization in the Euclidean case is done by introducing a set 
of local covariant boundary counterterms. These counterterms are needed not only for 
finiteness of the on-shell action [33^ [3^ [35l [22l [36] but also for the variational problem for 
AdS gravity to be well-posed [37j. In the Lorentzian setup, in addition to the infinities due to 
the non-compactness of the radial direction, there are also new infinities because of the non- 
compactness of the time direction. Correspondingly, in checking the variational problem one 
now has to deal both with boundary terms at spatial infinity and also at timelike infinity. 
Thus, in generalizing the Euclidean analysis to the Lorentzian case there are two issues to 
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be discussed. First, one has to check that the EucUdean analysis that leads to the radial 
counterterms goes through when we move from Euclidean to Lorentzian signature. This is 
indeed the case because all steps involved in the derivation of the radial counterterms are 
algebraic and hold irrespectively of the signature of spacetime. The second issue one needs 
to analyze are the infinities due to the non-compactness of the time directions and the new 
boundary terms at timelike infinity. A complete analysis of this issue requires knowledge of 
the asymptotic structure of the solutions near timelike infinity, which as far as we know is 
not available. Our prescription bypasses this problem by gluing in Euclidean AdS manifolds 
near timelike infinity. This effectively pushes the asymptotic region to the (radial) boundary 
of the Euclidean AdS manifold, whose asymptotic structure is well known. 

What remains to analyze is whether there are any problems at the 'corners', i.e. at the 
hypersurfaces where the Lorentzian and Euclidean solutions are joined. In principle, there 
can be new corner divergences which would require new counterterms. In this section we 
show that such corner divergences are absent in two examples: a free massive scalar field 
in a fixed background and pure gravity. We expect such corner divergences to be absent in 
general. 



This subsection serves to illustrate the problems at hand, and we will therefore adopt 
the simplest possible setting. As indicated in Fig. [21 we consider a single corner where 
the contour makes a right angle, passing from a vertical segment to a horizontal segment. 
The spacelike manifold is taken to be M'^^^. In the absence of sources, we explain below 
how this contour can be 'filled' with empty AdS^+i, with a metric that jumps along a 
spacelike hypersurface from Euclidean to Lorentzian. On this background, we consider a 
massive scalar field which propagates freely and without backreaction, and we compute the 
renormalized one-point function of the dual operator. 

3.1.1 Background manifold 

For the bulk manifold under consideration, we take one copy Mi of empty Lorentzian 
AdSd+i in the Poincare coordinate system (r, x*) with the metric 



3.1 Scalar field 



ds^ = dr^ -|- e^^riijdx 



(3.1.1) 



and one copy Mq of empty Euclidean AdS in similar coordinates and metric 



ds 



,2 



dr'^ + e^'^dijdx 



■'dx^ . 



(3.1.2) 
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drMo 

drMi 



t 

Figure 2: A single corner in the contour in the complex time plane. We use this part of a 
field theory contour to illustrate the holographic renormalization. 

We will take x*^ to be the time coordinate, denoting it by t on Mi and r on Mq. We use 
the notation for the other boundary coordinates, so for example = (t,x'*) on Mi, and 
we also introduce = (r, x"^). The conformal boundaries of the spacetimes lie at r — > cx) 
and are denoted drM\ and drMo. 

Next, we perform the gluing and the matching. To this end, we cut off the spacetimes 
across the surface t = and r = such that t > and r < 0, and glue them together 
along the cut surface which we call dtM. This surface is the extension of the corner in the 
boundary to the bulk. The induced metric on dtM is the same on both sides, 

hABdx^dx^ = dr'^ + e^"" 5abdx''dx^ , (3.1.3) 

and the extrinsic curvature ICab vanishes on both sides. Therefore, both the conjugate 
momentum ttab = Hab — K-ab as well as the induced metric Hab are continuous across 
dtM and all the matching conditions of section [2.2.11 are satisfied for this background. (We 
elaborate on the matching conditions for gravity in the next subsection.) The unit normals 
to dtM on either side are given by 

n[i]fidx^ = —e^dt , n^^dx^ = e^dr , (3.1.4) 

where we used subscripts in square brackets to indicate whether we are on Mi or on Mq. 
We will use this notation throughout the paper. 

Notice that the contour of Fig. [2] is not complete, since there is no out state specified at 
the right end of the contour. This should be remedied, for example by gluing a Euclidean 
segment at t = T which would result in the vacuum-to-vacuum contour of Fig. [T^. In the 
bulk, this incompleteness means that we should also glue another solution to some 'final' 
hypersurface lying in Mi. To obtain the contour of Fig. [1^, for example, one should glue in 
half a Euclidean solution M2. With M^"^ as the spacelike manifold, this would result in a 
spacetime as sketched in Fig. [7^. 
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In this section we will focus on a single corner. We will therefore omit any contributions 
from such an M2, as well as some terms defined on the final matching surface for Mi. Since 
the matching between Mi and M2 is a word-for-word repetition of the matching between 
Mq and Ml, these terms can be easily reinstated. 

3.1.2 Scalar field setup 

In the background we just described, we consider a scalar field $ of mass m, dual to a scalar 
operator O of dimension A such that = A{A — d). We will consider the case where 
A = d/2 + k with k G {0, 1, 2, 3, . . .}, and sometimes we will specialize to A; = 2. The actions 
for <I> on Ml and Mq are respectively given by: 



Suppose $ is a solution on Mq and Mi of the equations of motion derived from these actions, 
with asymptotic value corresponding to the radial boundary data and furthermore satisfies 
the aforementioned matching conditions (which we discuss in more detail below) on the 
gluing surface. Our aim is then to compute the corresponding on-shell action, 



while using the method of holographic renormalization to make it finite. Note that (I3.1.6P 
can alternatively be written as: 



with a path C in the complex time plane as in Fig. [2l which goes down at first (yielding 
—Sq after substituting t = —it) and then makes a corner and lies along the real t axis. We 
will not use this notation in this example, but it will be relevant when we consider gravity 
below. 

The holographic renormalization relies on the fact that the solution can (both on Mq 
and on Mi) be written as a Feffer man- Graham expansion: 



$ = e('=-'^/2)^-('/'(0) + e-'>(2) + . . . + e-2'^q0(2fc) + 0(2fc) loge-^T + ...). (3.1.8) 




(3.1.5) 



(3.1.6) 




(3.1.7) 





(3.1.9) 
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with □ the Laplacian of the boundary metric on d^M^ which in the case at hand is either 
r]ij or 5ij. Similarly, all coefficients 4'{2n) ^^4>{q) for n < k and 0(2fc) oc C\^(f)(^Q^, all with 
some fc-dependent coefficients. The coefficient 4'{2k) is normally nonlocally determined by 
(/)(o), but in our case it also depends on the initial data that one may specify at dfM. In 
Euclidean backgrounds without corners, this coefficient (times a factor —2k) is precisely 
the renormalized one-point function [22]. Below, we show this is still the case in Lorentzian 
signature and in the presence of corners. 

3.1.3 Matching conditions 

Let us first discuss the matching conditions of section 12.2.11 in more detail. Consider two 
solutions <&i and on Mi and Mq that satisfy the given radial boundary data, but have 
arbitrary initial data. The first matching condition is continuity of <J> across dtM, that is: 

$o(r = 0, r, x") = ^>i(t = 0, r, x") (3.1.10) 

for all r and x'^. 

To derive the second matching condition, we compute the on-shell action for $o ai^d '^'i 
satisfying the equation of motion and the first condition (jS.l.lOh . This action is divergent 
and we regulate it by cutting off the radial integrals at some large but finite tq. We then 
consider the variation of the regulated version of the total action p.l.6p as we vary the 
initial data ^{t = 0, r, x") and obtain 

6{iSi -So)= [ Vhe~''{-idt'^i - 9^$o)<5$i , (3.1.11) 

J dtM 

where we used that S^i = 6^o by (jS.l.lOp . As explained in section O we then request 
that the total action is also at an extremum with respect to the initial data. The second 
matching condition thus becomes: 

idt<^i + dr<^>o = on dtM . (3.1.12) 

As we mentioned before (and as one may check easily using t = —it), this second matching 
condition can be read as C^-continuity in the complex time plane of ^ across the corner. In 
the remainder of this section, whenever we write (7"-continuity, we always mean continuity 
in the complex time plane. 

Now let us substitute the Fefferman-Graham expansion (j3.1.8p of $i and <I>2 in the 
matching conditions (j3.1.10|) and (j3.1.12p . The matching conditions imply the C^-continuity 
of all coefficients 0(2/), which, in turn, implies higher-order continuity of the source (^(q)- 
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For example, the first matching condition for (j)(^2) becomes, via (|3.1.9|) . 

O[i](t>{o)[i] = 0[o]<P{om on dtM, (3.1.13) 

which shows that (j)^^-^ has to be at least C^-continuous across the matching surface. Notice 
that this is again continuity in the complex time plane, since □[!] is not equal to Dfo]- Next, 
the second matching condition applied to (j)^2) actually implies C^-continuity for 0(o): 

idtO[-^](f>(o)[i] + <9rn[o]</'{o)[o] = . (3.1.14) 
A similar story holds for the subsequent terms. Since the highest number of derivatives 
is always in 4'{2k) ^^4'{o): applying the second matching condition to this term results 
eventually in a C^'''''~^-continuity condition for 0(0) in the complex time plane. Below, we 
will see the relevance of these high-order continuity conditions. 

A comment considering this smoothness condition for (/)(q) is in order. Namely, this con- 
tinuity condition essentially follows from the requirement of the existence of a Fefferman- 
Graham expansion at the matching surface. In that light, this higher-order smoothness 
condition for 0(o) is not surprising, since without it the Fefferman- Graham expansion would 
fail even in the case without a corner. Although it would be interesting to study what 
happens for discontinuous boundary data, such an investigation can be undertaken inde- 
pendently of the presence of corners and shall not be pursued here. 

3.1.4 Holographic renormalization 

The on-shell action ()3.1.6p . evaluated on the solution that satisfies the matching conditions, 
is of the form: 



iSi- So = -^ [ ^/^^l^r^l j ^<^Qdr^Q 

^ Jdrhh ^ JdrMo 

-If Vh[-i<i>idt<i>i + $o5r^o] ■ 

^ JdtM 



(3.1.15) 



The contributions from dfM, i. e. the second line in (|3.1.15|) . vanish by virtue of the matching 
conditions. Recall that we are omitting the contribution from any 'final' surface for Mi, 
which will however by the same mechanism cancel against a matching solution. 

The remainder of the action is defined on the cutoff surface r = ro and it would diverge 
if ro ^ cxD. Therefore, a counterterm action has to be added before removing the cutoff. 
Since the radial terms in (j3.1.15p have a familiar form, one can use the usual procedures 
of holographic renormalization to find the counterterm action [6]. Let us for example take 
k = 2, for which 

= 5 L ("= - |)*^ + ^) + ^«?*.og.-) (3.1.10) 
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is the counterterm action. The first two terms are actually valid for any k > 2 and we used 
the notation for the Laplacian of the induced metric 7 at r = tq. In our case, we simply 
have = e~^''n, both on Mi and on Mq. Taking care of the signs, we find that 

iSi - So + iSct,i + Sctfi (3.1.17) 

is finite as rg — > 00. We see that the usual procedure of holographic renormalization yields 
a finite on-shell action and possible initial or final terms (which might have caused corner 
divergences) are absent exactly by the matching conditions. 

3.1.5 One-point functions 

One-point functions are computed by taking variational derivatives of the on-shell action 
with respect to the boundary data. Let us compute the one-point function {0[ij{x)), where 
the subscript indicates that x lies on drMi. In QFT on a background with a Lorentzian 
metric g(o)ij, the coupling between a source (/)(o) and an operator O in the partition function 
is as in (12.1.ip . Therefore, the one-point function is 

(0[i](x)) = ^^^_Z[</.(o)]. (3.1.18) 

In our case, the partition function Z[(^(o)] is given by the renormalized on-shell supergravity 
action. The easiest way to take care of the divergences is by taking the functional derivative 
before removing the regulator, resulting in: 

• (fc+d/2)ro X p . 

(0[i](x)) = hm iSi -So + iSct,i + Sct,o , (3.1.19) 



where the extra factor ^C^^+d-f^yo converts $ to </)(o) and 7 to (^(o) as ro ^ 00. 

In performing this computation, we see that the presence of corners gives rise to corner 
terms, which arise from the integration by parts that is necessary in varying the counterterm 
action (j3.1.16p . For example, for the variation of the second term in (|3.1.16p we obtain: 

1 f ^ <^\J^^ \ _ f /rT'^^°7^ , 1 f /nre-2'-(5t$5$ - $9(5$) 



.2 7a^M ^ "'2(1-A:)y y^^jv/ ^ ' ' 2(1 - A:) 2 ^ ' ' 2(1 - k) 

The second term on the right hand side is a corner contribution. However, a similar corner 
term arises in S'ct,0) and in the total action ()3.1.17p these two corner terms cancel each other 
precisely by the matching conditions. 

The subsequent terms in the counterterm action are all of the form ^/7$n"$ for n < k, 
plus a log term of the form ^/7$n^$ log e~^°. After the integration by parts, these all give 
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corner terms as well, which involve a higher number of derivatives of More precisely, the 
corner expressions that one obtains from such terms are of the form 

V7e~'"^5tn"-i$ , (3.1.20) 

c 

and equivalent terms with some of the derivatives shifted to the first 

Let us now systematically show that all such terms cancel against a matching solution, 
using the higher-order smoothness of (^(q) that we derived before. First of all, recall that 
the matching conditions imply that 0(o) should actually be C^'^'^^-continuous. This in turn 
means that (^(2) is C^*^~^ continuous, ^(4) is C^*^~^-continuous, etc., up to (t)[2k) and (/)(2fc), 
which are just C^-continuous. Substituting this in the Fefferman-Graham expansion p.l.Sp . 
we see that <I> is not only C^-continuous by the matching conditions, but also C'^-continuous 
up to terms of order ^-i^+'^/'^y ^ and continuous up to terms of order g(-'=-'^/2+2)r^ ^^^^ 
We now rewrite the leading piece of ()3.1.20p as 

I ^(,+,/2-2n)r (a,□^^$) + . . . (3.1.21) 

Jc 

A complete cancellation of this term between Mi and Mq takes place if <l> is C^"~^- 
continuous up to and including terms of order e~('^+'^/^~^")''. However, the previous ar- 
gument shows that C^"~^-continuity for $ holds up to terms of order g-('=+<^/2-2n+4)r^ 
the continuity condition is satisfied indeed, for all n < k. Therefore, as tq — > oo, the terms 
coming from Mq and Mi cancel indeed and no corner contributions to the one-point func- 
tions arise. A similar argument shows that there is no problem with the log term with 
n = k either. 

Having shown the absence of corner contributions in ()3.1.19p , one finds that the expres- 
sion for the one-point function becomes of the standard form, given for example for k = 2 
by: 



ro-»oo 



r=ro 

(3.1.22) 



Substitution of the expansion (j3.1.8p yields the familiar result: 

(0[i](x)) = -2/c0(2fe)[i](x), (3.1.23) 

which is actually valid for all nonzero k, see for example [6]. 

Finally, consider the one-point function on M2, where we should use the Euclidean 
version of the source-operator coupling, — / y^^Q^(/)(0)O- Repeating the above procedure, 
we find again: 

(C[o](a;)) = -2/c(/>(2fe)[o](x). (3.1.24) 
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Since 4'(^2k)i^) is continuous across the matching surface by the first matching condition, 
and since locahzed corner terms are absent, the one-point function is continuous across the 
corner as well. 

3.2 Gravity 

For gravity the procedure requires modification and becomes more involved. We therefore 
begin with an outline of the steps taken below. 

The first step in this procedure is to establish the variational principle for the Einstein- 
Hilbert action for a manifold whose boundary has corners. Recall that in the Euclidean 
setup a well-defined variational problem requires the addition of the boundary counterterms 
|37j and the variational derivatives w.r.t. boundary data lead to the boundary correlators. 
In the Lorentzian setup the variational derivatives w.r.t. initial and final data are also 
important and lead to matching conditions. The analysis of the variational problem is done 
in subsection I3.2.2[ We will find that there is a need for a special corner term. 

The next step is to understand how to glue the various pieces together. Given a cor- 
ner in the boundary contour there should exist a corresponding bulk hypersurface across 
which the various bulk pieces are matched. So we need to understand the possible bulk 
extensions of the boundary contour. This is analyzed in subsection 13.2.31 where we show 
that the extensions are parametrized by a single function /(r, x") with a certain asymptotic 
expansion. 

Using these results we then derive the matching conditions in subsection 13.2.41 and find 
their implications for the radial expansion of the bulk fields near the corner in 13.2.51 These 
are all the data we need to analyze whether there are any new contributions to the on- 
shell action and the one-point function from the matching surfaces. This is done for the 
on-shell action in subsection l3.2.6l and for the 1-point functions in subsection 13.2.71 We find 
that there are possible contributions from each segment but the matching conditions imply 
complete cancellation between the contributions of the two pieces that one glues to each 
other. 

The upshot of the discussion is therefore very similar to the scalar field: we will show 
that no localized corner terms arise and that the one-point function of the stress energy 
tensor is (appropriately) continuous across the corner. 
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3.2.1 Setup 

As we mentioned earlier, we consider manifolds Mq consisting of a number of segments Mj 
where the metric is Lorentzian or Euclidean. To simplify the computation of the on-shell 
action for these spacetimes, we introduce a notation where the Einstein-Hilbert action Sj 
for each separate segment Mj is always written as 

Sj = ^ [ d'^+^x V^{R - 2A) , (3.2.1) 

where k? = 87rGd+i and A = —d{d — l)/(2^) with £ the AdS radius. Throughout this 
paper, we set i = 1. In p.2.ip the square root is defined with a branch cut just above the 
real axis. For example, for a Euclidean metric \J —G = —i\/\G\ , so that iS = —Se with 
Se = J \/\G\ {—R + 2A) the correct Euclidean action. Similarly, for a Lorentzian metric on 
a backward-going contour we obtain an extra minus sign since we are on the other branch of 
the square root. (To see this, notice that the time coordinate tc on this segment is given by 
t = e^'^tc- If G, Gc denote the metric determinant in the t, tc coordinate system, respectively, 
then Gc = e^'^*G, and we make a full turn indeed.) The advantage of this formalism is that 
the total Einstein-Hilbert action Seh for Mc becomes 

iSEH = iSo + iSi + ... (3.2.2) 

for all vertical or horizontal segments Mq,Mi, . . . We see that all the signs are absorbed 
in the volume element. This action for Mq needs to be supplemented with various surface 
terms which we define in due course. 

Although we will not discuss this in detail, this prescription can be extended to general 
complex metrics, allowing for the 'filling' of more general QFT contours that are not just 
built up from horizontal and vertical segments in the complex time plane. In such cases 
the bulk metric G^u may be complex, but it should always be non-degenerate for the 
scalar curvature to be well-defined. Allowing for a complex metric implies that one has to 
allow for complex diffeomorphisms as well, for example to bring the metric to a Fefferman- 
Graham form. Complex diffeomorphisms are discussed in some detail in [32]. For such 
cases, our choice for the branch cut in the volume element is then precisely consistent with 
the requirement that a QFT contour cannot go upward in the complex time plane. 

3.2.2 Finite boundaries 

In equation (j3.2.2p . we split the on-shell action for Mq as a sum over the various segments 
Mj. Just as for the scalar field, we will find the matching conditions via a saddle-point ap- 
proximation which involves taking functional derivatives of the on-shell action with respect 
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to the initial and final data. This only works if we have a well-defined variational principle 
for each segment separately, which is what we investigate in this subsection. 

Consider a single Asymptotically locally AdS (AlAdS) manifold M with a (possibly 
complex) metric G^i, and two 'initial' and 'final' boundaries which we denote here as d±M. 
The manifold M also has a radial conformal boundary, which we denote as drM, and the 
corners where d±M meets drM are denoted as C±. We pick coordinates (r, x*) on M, with 
= {t, rr"), and we will also use x"^ = (r, x""). The conformal boundary is again at r ^ oo. 
We regulate the computation of the on-shell action by imposing r < rg. In this subsection 
we consider the variational principle in the case where one keeps tq finite throughout. 

A well-defined variational principle for Dirichlet boundary conditions in the presence 
of corners requires the Einstein-Hilbert action to be supplemented not only with the usual 
Gibbons-Hawking- York boundary terms on drM and d±M, but also with special corner 
terms defined on C± [38l ESI HD] • To find these corner terms, we choose coordinates such 
that drM is given by r = tq and d±M hy t = t^-. The metric near the corners can be put 
in the following two ADM-forms: 

G^.ydxi'dx'' = {H^ + HiH')dr^ + 2Hidx^dr + %dx'dx^ , 
%dx'dx^ = (-M^ + MaM^)dt'^ + IMadx^'dt + aabdx^'dx^ , 

as well as 

G^^dx^dx" = (-M2 + MAM^)df + 2MAdx^dt + hAsdx^dx^ , 
hABdx^dx^ = {H"^ + HaH'')dr'^ + IHadx^dr + aabdx^'dx^. 
Relating the two metrics, we find 



(3.2.3) 



(3.2.4) 



= ^ Ht = Mr (3.2.5) 

M^ = M''-{M^-fH^ -ML = ^, 

For a real Lorentzian metric and are positive, whereas they are negative for a 
Euclidean metric. We will henceforth assume that o", the determinant of fJafe, is real and 
positive. This will simplify the discussion and is sufficient for all the examples below. 

The standard Gibbons-Hawking- York surface terms involve the extrinsic curvature ^KLab 
of d±M and Kij of drM, which we will define using the (possibly complex) unit normals. 



drM : hf.dx'' = ^^^=dr Kij= — ^= — {DiHj + DjHi - dr%) , (3.2.6) 

v-7 2^J^pm 



d±M: ^n^dx'' = ±^^dt ^ ±}CAB=±^^^^-{DAMB + DBMA-dthAB). 
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Adding the Gibbons-Hawking- York terms, we define the bare action as: 
1 



j (f^^x^/^{R-2K) + 2 f + 2 [ d'^x^/^K], (3.2.7) 

J Jd±M JdrM J 

where here and below the summation over d±M is imphcit and we use the conventions 

of subsection 13.2.11 for the square roots. For a real Lorentzian metric all the above terms 

are real, but for a real Euclidean metric all terms in (I3.2.7P are purely imaginary (because 

V— G and are then imaginary, from which it follows that ^n^^dx^ and therefore ^K, are 

imaginary as well) . As one may verify explicitly, in the latter case our choice of branch cut 

for the square roots in the volume elements implies that iSh = —Se with Se the Euclidean 

action with the correct Gibbons-Hawking- York terms. 

In the case of corners, ()3.2.7I) is not the correct action to use for Dirichlet boundary 

conditions. This is because we cannot perform a diffeomorphism at the corner mixing t and 

r without changing the definition of the two slices and therefore Ht, Mr, and are no 

longer pure gauge at the corner. With this in mind, the variation of the bare action (I3.2.7P 

is given by the equations of motion, the conjugate momenta for all the various boundaries, 

plus a corner term 



2^2 

with X given implicitly by 



§Sh = ^ [ d'^-^xy/a5X± + . . . , (3.2.8) 



SX± = ±2 . (3.2.9) 

To find an explicit form of X±, we have to integrate 6X for fixed ^nd H^, using the 
relations p.2.5p . If the metric is completely real and and are positive, then we find 

5X± = ±2(5arcsinh(^^^j , (3.2.10) 

whereas if negative and and are positive we get 

(HAP \ 

) . (3.2.11) 

We can rewrite these expressions in a covariant form using the unit normals defined in 
(|3.2.6p . Their inner product is given by: 

^n^'n^ = i-^iW" . (3.2.12) 



For real M^, ^/IP and M^, we can therefore write without branch cut ambiguities: 

2 arcsinh( ^n^^h^) JVP > 



(3.2.13) 

2i arcsin(i='=n^n^) Af^ < Q. 
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In the more general case, the required corner term has the same structure but one needs to 
be careful about the branch cuts. Notice that X is defined up to a local piece, for example 
a constant. 

Following [381 [39l BO] , we aim for a variational principle that is well-defined for a fixed 
induced metric on the boundaries, i.e. for fixed 7jj and Kab- In that case, we should add a 
corner term to cancel the unwanted variation 5X in ()3.2.8p . Such a corner term is given by 



(3.2.14) 



Adding corner terms to the action ()3.2.7p defines an improved (but still bare) action 5/, 



Si 



Sb + Sc^ 
1 



2^2 



(f+^xy/^{R - 2A) + 2 / d^x^/h^K 

Jd±M 

+ 2 [ d'^Xy^k- I d'^-^X^/^X± 
JdrM Jc± 



(3.2.15) 



whose variation is of the form 



^JdrM Jd±M 



d'^~^x6{^)X± 



c± 



(3.2.16) 



which is the correct variation for Dirichlet boundary conditions indeed. We will henceforth 
use this improved action as the bare action and drop the subscript I. 

3.2.3 Fefferman- Graham coordinates 

The above discussion was valid for a general spacetime whose boundary has corners. Since 
we are interested in AlAdS spacetimes where the metric diverges near the radial boundary, 
we will run into divergences as we let rg oo. To investigate these divergences, we pick a 
coordinate system in which the metric is of the Fefferman- Graham form. 



ds'^ = dr"^ + 'jijdx''dx^ 



with the radial expansion 



-2r 



9{2)ij 



+ e '^''[9{d)ij + 9(d)ij log e + • 



(3.2.17) 



(3.2.18) 



From the Einstein equations we find that all coefficients g(2n)ij with 2n < d, as well as 
9{d)ijj locally determined by g[o)ij, and involve up to 2n or d derivatives of <7(o)ij- The 
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term g(^d)ij is not locally determined (except for its trace and its divergence) and this term 
directly enters in the one-point function of the stress energy tensor |22j . 

The disadvantage of the Fefferman-Graham form of the metric is that one can generally 
no longer pick a coordinate t such that the surfaces d±M are given by slices of constant 
t. On the other hand, one can use the leftover gauge freedom to make sure that d±M are 
asymptotically given by: 

5±M:t = /±(r,x'^), (3.2.19) 

with 

lim /±(r,x'') = t± (3.2.20) 

and t± constants. We will discuss the asymptotic behavior of f± more precisely below. 

Let us consider a single initial or final boundary. Dropping for now the subscript it, we 
write an ADM-decomposition of ^ij near the corner: 

-iijdx'dx^ = {-N^ + N''Na)dt'^ + 2Nadtdx°- + Tabdx^dx^ . (3.2.21) 

We may pick boundary Gaussian normal coordinates centered at the corner, so that Na ~ 
0(1). Furthermore, iV^ = e'^^'N'^^-^ + N^^) + ■■■ and Tab = e^'"r(o)afe + T(^2)ab + ... We can 
relate this ADM-decomposition to the double ADM-decomposition (j3.2.4l) of the previous 
subsection by introducing a new coordinate 

t' = t-fir,x''), (3.2.22) 

after which the initial slice is given by t' = 0. In the new coordinates, the metric is of the 
form (j3.2.4p . with t replaced by t' , and with the components 

-M^ + MaM^ = -N^ + NaN"" 

Mr = {-N^ + iVaiV")9r/ 

Ma = Na + {-N^ + N,N^)daf 

(3.2.23) 

+ HaH" = 1 + i-N^ + NaN''){drff 

Ha = {-N^ + N,N')dafdrf + Nadrf 

(Jab = Tab + {-N^ + N,N'')dafdbf + NaObf + NbOaf , 

where indices are raised with the appropriate metric. We use these equations below to 
write down a radial expansion of the components on the left-hand side in terms of the 
Fefferman-Graham expansion and a radial expansion of /. 

For AlAdS spacetimes the Dirichlet boundary data are given by 5(o)jj ™d h^s- Asymp- 
totically, g(o)ij determines a Fefferman-Graham radial coordinate as well as the subleading 
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coefficients up to g(d)ij the FefFerman-Graham expansion of the metric. Of course, the 
initial and final metric Hab should be such that d±M can be embedded in the asymptotic 
metric dictated by 5(o)ij and this condition constrains the asymptotic form of /iab- To be 
precise, h^s should have a radial expansion that is compatible with the last three equations 
in (I3.2.23P for a certain /. However, as long as / is unspecified, Hab is not to any order 
determined in terms of g(o)ij ■ 

We remark that the last three lines in (j3.2.23p signify constraints on h^B only. Therefore, 
they are different from the usual constraints on the initial data in a Hamiltonian formalism 
of general relativity, which also involve the extrinsic curvature. These usual constraints are 
satisfied if the extrinsic curvature of the initial slice is computed using the embedding of the 
initial slice as a hypersurface in the solution. Therefore, they are automatically satisfied if 
we compute the extrinsic curvature using the first three lines of ()3.2.23p . Since this is how 
we compute the extrinsic curvature below, we will not worry about these constraints. 

3.2.4 Gluing and matching conditions 

In the previous subsections, we found an improved action (j3.2.15p and discussed the Fefferman- 
Graham expansion for a single AlAdS spacetime with corners. We now take two of such 
spacetimes and glue them together along the initial and final hypersurfaces d±M. 

We will denote the two segments by Mq and Mi and we glue S+Mq to d-Mi, which we 
from now on we denote as dfM. The corner, i.e. the intersection of 9jM with drMo and 
dr-Mi, is denoted by C. As before, a subscript (sometimes in square brackets) indicates the 
manifold under consideration. We make no assumptions about the signature of the metric 
on Mq , All and in fact the metric may even be complex. We write the total action as 

iSo + iSi , (3.2.24) 

with the individual actions given by ()3.2.15p . We recall that we use the conventions of 
subsection l3.2.H so extra factors of i might be included in the volume elements and extrinsic 
curvatures. As we did for the scalar field, we will henceforth ignore the contribution from 
other segments than Mq and Mi as well as the contribution to the on-shell actions of Mq 
and Ml that may arise from other matching surfaces. 

Let us now find the precise matching conditions that the metrics on Mq and Mi have to 
satisfy near dtM. The first matching condition is continuity of the initial and final Dirichlet 
data. For gravity, this becomes continuity of the induced metric: 

h[Q]AB = h[i]AB ■ (3.2.25) 
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The second matching condition is obtained from the variation of the on-shell regularized 
action with respect to the data on dfM. Let us first suppose the variation vanishes at 
the corner C. In that case, we read off from (|3.2.16|) that the second matching condition 
becomes: 

^[o]AB + = . (3.2.26) 

We can also consider a variation that does not vanish at C, for which (13.2. 16p shows that 

(X[i]+X[o])5(^/^) = 0, (3.2.27) 

where we included the S{^/a) because of the following reason. Notice that this is a corner 
matching condition which is therefore not valid to all orders in r. However, since ^/o^ ~ 
g{d-i)r^ (|3_2.27p is actually divergent as tq — > oo. Therefore, it only vanishes completely if 
the X's match to high order in their radial expansion. If there are no log terms, then we 
find 

+ X[o] = 0(e-'^'-) . (3.2.28) 

Equation (I3.2.28P is the corner analogue of the bulk matching condition (j3.2.26|) . Notice 
that such a corner condition was absent when we discussed the scalar field discussed above. 
Its implications will be investigated in the next subsection. 

We showed before that ICab and X are imaginary for a Euclidean metric. Therefore, 
although it is not transparent in our notation, the matching conditions (j3.2.26p and (I3.2.28P 
do contain factors of i when joining a Lorentzian to a Euclidean metric. 

3.2.5 Imposing the matching conditions 

For the scalar field, the matching conditions were crucial in demonstrating the cancellation 
of corner divergences and the absence of localized corner contributions to the one-point 
function. A similar cancellation will occur for gravity, but imposing the three matching 
conditions ()3.2.25p . ()3.2.26p and ()3.2.28p will not be as straightforward as for the scalar 
field. 

In this subsection we shall impose the matching conditions order by order in a radial 
expansion of /iab, K^ab and X. We start with a detailed analysis of the leading-order terms 
in the matching conditions. We then discuss continuity in the complex time plane of the 
boundary metric. Just as for the scalar field, the higher-order continuity is related to the 
continuity of the subleading terms in the Feffer man- Graham expansion of the bulk metric. 
Afterwards, we show that our leading-order results extend to the higher-order terms as well. 
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Leading order matching conditions 

We will work in the FefFerman-Graham coordinates, with the matching surface dtM given 
by (|3.2.19|) . Without loss of generality, we assume that the corner is given by t = on dj-Mi 
and T = on drMo, so linir^oo fir, x"") = on either side. We suppose that / behaves 
asymptotically as 

/ = e-'^/(i)(x'^) + ... (3.2.29) 

This is the leading asymptotic behavior of /, since any slower falloff near r — > oo would 
yield a non-spacelike induced metric on dfM in a real Lorentzian spacetime. Substituting 
(j3.2.29p and the leading-order terms in the ADM-decomposition p.2.2ip of jij in (|3.2.23p , 
we find the leading behavior of H^, and M^. The inner product between the unit 
normals, given in (j3.2.12p . becomes to leading order: 



T Z'^'"^'" (3,2^30) 

Since continuity of X± follows from continuity of ^n^fi^, the corner matching condition 
(j3.2.28p becomes to leading order: 



\/^(o)[o]/(i)[o] = \l^{Q)[i]f(im ' (3.2.31) 

where we reinstated the subscripts to indicate the manifold under consideration. 

Let us work out the consequences of this condition. Recall that we absorbed factors 
of i in the square roots of the metric determinant, and therefore ()3.2.3ip is not necessarily 
a relation between real quantities. For example, if we match a Lorentzian to a Euclidean 
solution, then N'^q^ changes sign across the corner and the square root on the Euclidean side 
of (|3.2.3ip becomes imaginary. On the other hand, the square root on the Lorentzian side 
is real, and so is /(i) since we use real coordinates. This means that in that case we must 
have: 

/(i)[o] = /{i)[i] = 0, (3.2.32) 

which more generally holds in all cases for which the phase of ^f^-^ is discontinuous across 
the corner. Actually, this phase is only continuous when we match two solutions with 
the same signature (recall that we chose boundary coordinates in which A'^a(o) = 0)- This 
happens either if we have no corner at all, or if the corner makes a 180-degree turn. In the 
first case, we can pick boundary coordinate systems in which A^^^^ is continuous across the 
corner and ()3.2.3ip becomes simply 

/(i)[o] = /{i)[i] • (3.2.33) 
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Figure 3: On the left, the dotted hnes represent two bulk hypersurfaces given by t = / in 
the vicinity of a corner in the boundary contour at t = 0. On the right, we see that around 



Since we just artificially split a spacetime in two parts, it is natural that there is no further 
constraint on /. The case in which the corner makes a full turn is slightly more involved. 
First of all, the two boundary segments ending on the corner must be straight horizontal 
lines in the complex time plane, since the boundary contour cannot go up in this plane. 
We may again assume that N'^^-^ is continuous across the corner, but that does not mean 
that the square roots in (|3.2.3ip are. Namely, one of the segments is backward-going in the 
complex time plane and in subsection 13.2.11 we already mentioned that the square root for 
a backward-going contour results in a minus sign. The matching condition for a full turn 
therefore becomes 



This implies that, at least at this order, we can freely move the hypersurface dtM up and 
down in the bulk, as long as we move it by the same amount on both components and keep 
the location of the corner fixed. We have sketched this in Fig. [3l 

We have worked out the leading order term in the corner matching condition in three 
cases, corresponding to three different corners. We emphasize that our formalism of subsec- 
tion 13.2.11 allowed us to summarize all three cases in the single equation (j3.2.3ip . We will 
see below that the subleading behavior of / is constrained in an analogous way. 

As a sidenote, let us also compute the leading order term in the radial expansion of the 
second matching condition (|3.2.26p . If we use (j3.2.6p to expand the trace of the extrinsic 
curvature ^JCab the leading order term becomes: 



a full turn in the boundary contour it is natural to expect that /[i] 



/[o]- 



/(1)[0] 



/(i)[i] ■ 



(3.2.34) 




(3.2.35) 
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The trace part of (j3.2.26|) therefore results to leading order again in (|3.2.31|) . It is plausible 
that for AlAdS spacetimes the corner matching condition ()3.2.28p follows from (I3.2.26P 
and does not need to be imposed separately. This would be related to the fact that the 
asymptotics of the bulk metric are completely determined by the Fefferman-Graham data, 
but a more complete analysis is required to settle this issue completely. This will not be 
attempted here and we will instead continue to treat ()3.2.28p as an additional condition. 

Continuity in the complex time plane 

Just as for the scalar field, the Fefferman-Graham expansion relates subleading terms in 
the matching conditions to higher-order continuity in the complex time plane of the sources. 
Before proceeding with the subleading terms in the matching conditions, let us therefore 
first discuss the notion of smoothness in the complex time plane for the boundary metric. 

Consider a contour in the complex time plane with a corner. We define C'^-smoothness 
for the boundary metric as the condition that the /c'th order t-derivatives of the metric 
components exist at the corner of the contour. Although this is a natural definition, in 
our notation a complication arises because we do not work directly with a complex time 
coordinate on for example the vertical segments. Instead, we rather use a contour time like 
tc or T which is real on a particular segment of the contour and for such parametrizations 
the continuity condition has a different form. We may find this new form by regarding 
these local parameters as related to t via a complex diffeomorphism, for example t = —it or 
t = 2T — tc- If we use these parameters to express continuity of the metric, then we need to 
take care of the transformation properties under the diffeomorphism as well. For example, 

continuity of g{o)ij across the corner of Fig. [21 where t = —ir, becomes the condition 
that at the corner 

5{0)[0]rr = -9{0)[l]tt , 5(0)[0]ra = -^5(0)[l]ia 5 ff(0)[0]afe = 9{0)[l]ab ■ (3.2.36) 

Similarly, continuity in the complex time plane becomes dr{g(o)[o]ab) = ~'''dt{g(o)[i]ab) 
and (?T-((7(o)[o]rT) = ^f^t(5(o) [!]«)• The extension to higher orders and other components is 
analogous. As an example, take ds^ = dr"^ + 5abdx°'dx^ and ds^j^j = —dt^ + 6abdx°'dx^ . 
Although there is an apparent discontinuity in the metric components, with our definitions 
the metric is C°° at the corner. 

We will from now on assume that the boundary metric at the corner is continuous in 
the complex time plane. The reason for this smoothness condition is the same as that for the 
scalar field: it guarantees the existence of a Fefferman-Graham expansion of the metric at 
the corner, and the locally determined coefficients in this expansion are then automatically 
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continuous across the corner as well. Since we continue to use real coordinates like r, we 
will always need to supplement the continuity condition with the transformation under the 
complex diffeomorphism. 

Higher order matching conditions 

We showed above that the leading order matching conditions imply that /(^^ usually 
vanishes, except in special cases when N'^^^s^ does not change across the corner. In this 
subsection, we show that the matching conditions and the continuity of the boundary 
metric fix the higher-order terms in / to behave just as /(i), at least up to terms that vanish 
faster than e~'^^ . 

We first assume that the leading order term in / is: 

/(r,rE") = e-"7(n)(x'^). (3.2.37) 

One may easily check that in this case 

^n^h^ = Tyiv|^/(n)e(^-")^ + . . . (3.2.38) 

and a repetition of the previous analysis shows that, for n < d, the leading order term in 
p.2.28p becomes equivalent to 

\/%mf{nM = \/ %)[i]f{n)[i] ■ (3.2.39) 

Therefore, if the phase of is discontinuous across the corner, we find that not only 

/(I) but all terms up to and including order e~'^^ in f{r,x"') vanish as well. 

If N^Q^ is continuous, then /(i) does not necessarily vanish, equation (13.2.37P no longer 
holds, and the above derivation for the subleading terms is no longer valid. However, 
the continuity of the boundary metric implies that the locally determined terms in the 
Feffer man- Graham expansion ()3.2.18p are continuous across the matching surface as well and 
the metric is thus the same to high order on either side (up to the complex diffeomorphism 
discussed above). A discontinuity in ()3. 2.180 may appear at the earliest for the nonlocally 
determined term g[d)ij-, which is at overall order e^^"*^^*" in the radial expansion of the bulk 
metric. By substituting this radial expansion in the fourth equation in (I3.2.23p . and using 
the continuity to all orders of H'^ + HaH^, we find that / has to be continuous across the 
corner up to and including terms of order e~'^^. (Notice that the fourth equation in (j3.2.23p 
is invariant under / <^ — /, but we fixed the overall sign already at leading order.) 

This finishes our discussion about imposing the matching conditions: the previous two 
paragraphs show that / 'matches' up to and including terms of order e~'^^ for all three 
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cases. Up to this order, we find that /[q] = ~f[i] for a full turn, that /[qj and /[^j both 
vanish for any other corner, and that /[q] = /[i] if there is no corner at all. In the next 
subsection, we will use these conditions to demonstrate the absence of localized (divergent) 
corner contributions to the on-shell action, in order to eventually show the continuity of the 
one-point function of the stress energy tensor around the corner. 

3.2.6 Computation of the on-shell action 

The bare on-shell action (13.2. ISp has the usual Gibbons-Hawking- York contribution from 
dtM as well as an extra corner contribution. However, the matching conditions directly im- 
ply that these terms cancel between the two spacetimes. The total action (j3.2.24p becomes: 

iSo + iSi = At / d'^+^xy/^iR -2A) + A^ I d'^x^K 
J Mo JdrMo 

+ / d'^+^xV^{R-2A) + A^ [ d'^xy/^K. (3.2.40) 



k2 Jq^j^j^ 

This action can again be renormalized with the usual radial counterterms, except for a 
subtlety involving the bulk integrals in this action. Namely, the t-integrals do not run 
between fixed endpoints, say t = and t = T, but now rather end on t = f±{r,x°'). The 
usual radial counterterms, however, assume an r-independent limit on the bulk integral and 
the radial counterterms may not exactly cancel all divergences. 

We will now show that these extra divergences also cancel between the two matching 
solutions. To first order, the cancellation can be shown very explicitly. Namely, if / is of 
the form ()3.2.29p . then we can radially expand the volume element as: 

Gd'^+^x= [ ° dr [ dx" / dt^/N^ (3.2.41) 

Mo J J Jf(r,x^) 



ro 

dr 



dx" dtVN^ - e'-'^"^^' / dx"/(i)A/A^(o)^(o) 



+ 



The first term has an r-independent lower limit on the t-integral and so all divergences in this 
term are dealt with by integrating the usual radial counterterms also until /(ro,x"). The 
second term is not cancelled by counterterms and may lead to new divergences. However, in 
()3.2.40p a similar term comes from the expansion of the action Si for Mi and by the corner 
matching condition (I3.2.3ip the terms exactly cancel each other. Notice that an extra sign 
on Ml arises because we expand the upper rather than the lower limit of the t-integral. 

For higher orders, we recall that / is continuous or vanishing up to and including terms 
of order e"'^''. Using also the higher-order continuity of the bulk metric, a continuation of 
the expansion (|3.2.4ip shows that the corrections cancel up to finite terms. This means that 
no extra divergences arise from the discrepancy between the limits of the t-integration. 
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Having eliminated all possible sources of corner divergences, we may conclude that the 
usual radial counterterms are sufficient to make the total on-shell action finite. For example, 
in d = 4 the counterterm action is of the form: 



where the curvatures are those of the boundary metric 'jij at r = tq. This counterterm 
action is valid for all signatures if we define ^—7 in the same way as \/—G above, i.e. with 
the branch cut above the positive real axis. 

3.2.7 Continuity of the one-point function 

We have shown that the on-shell action can be holographically renormalized with the usual 
counterterms in the presence of corners. It remains to show that the one-point function is 
appropriately continuous around the corners as well. 

The renormalized one-point function of the stress energy tensor is obtained by varying 
the renormalized on-shell action with respect to radial boundary data. As for the scalar 
field, the integration by parts in the variation of a counterterm action like (|3.2.42p may result 
in localized corner contributions to the one-point function. However, a similar analysis as 
for the scalar field shows that the higher-order continuity of the boundary metric in the 
complex time plane ensures that such contributions again cancel between two matching 
solutions. 

Let us explicitly show the cancellation of the first corner term that arises from the 
integration by parts of the radial counterterms, which originates from the second term 
in ()3.2.42p . This is just an Einstein-Hilbert like term and it cancels against the matching 
solution if the extrinsic curvature of the corner, which we denote ir(o)afe) is continuous across 



Cancellation of the next term gives a higher-order continuity condition. Explicitly, the 
variation of these terms gives 




(3.2.42) 



the corner: 



K{0)[0]ab + -f^(0)[l]afe = . 



(3.2.43) 




+ / d''~'x^\niP'\V5jij - 7'Vj57fc/) + {V,P'^){n,j'''6jki - n'^Sjkj)] , (3.2.44) 



Jc 



where 



(3.2.45) 



4((i-l)(d-2) ((i-2)2 
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and n* is an appropriately defined unit normal for the corner as a submanifold of drM. 
From ()3.2.44p we explicitly see that the higher-order continuity condition involves up to 
three derivatives of the metric in d = 4. 

By the absence of initial or corner contributions, the holographic expression for the one- 
point function of the stress-energy tensor is completely analogous to the Euclidean case. In 
particular, it is expressed directly in terms of g(d)ij ^^id terms that are determined locally 
by g(o)ij- For example, in d = 4 we find up to scheme-dependent terms that 

i^ij) = ^{am - ^I(Trff(2))^ - T>5^2)] - Ii9f2)h + \9{2)ijT^ 9(2)) , (3.2.46) 

see |22] for the exact expressions in other dimensions. Alternatively, we may use the 'radial 
Hamiltonian' approach to holographic renormalization [361 which provides a more ef- 
ficient way of obtaining renormalized correlators. In this approach, the one-point function 
can be more compactly written as 

{T^J) = T^(d)ij , (3.2.47) 

where T^[d)ij is the term of dilatation weight d in the expansion of the radial canonical 
momentum in eigenfunctions of the dilatation operator. 

Since by assumption all locally determined terms in the Feffer man- Graham expansion of 
the metric are continuous, continuity of the one-point function will follow from continuity of 
9[d)ij across the corner. Fortunately, the continuity of g(d)ij follows directly if we substitute 
the expansion (j3.2.18p in the last equation of (|3.2.23p . The left-hand side in this equation 
is continuous to all orders by the first matching condition. On the other side, we know 
that / is continuous up to and including terms of order e"*^^, and we know that all g{2n)ij 
with 2n < d as well as g[d)ij are continuous since they are locally determined by g(o)ij- 
Collecting terms of overall order e^'^~'^^^ then establishes that g(d)ij has to be continuous as 
well. (As shown in |37j . there is no diffeomorphism freedom at this order if we fix a boundary 
coordinate system and a boundary metric, so continuity of g(d)ab implies continuity of g(d)ij 
indeed.) We have thus established that the vev of the stress-energy tensor is continuous 
across the corner (in the sense discussed in subsection 13 . 2 . Sh . 

We end this section with a remark about the function /(r, a;*^). Recall that we could 
in some cases freely specify this function at the corner, provided it was the same on both 
sides (possibly up to a sign). On the other hand, this function has no place in the QFT, 
and therefore holographically computed QFT correlators should be independent of /. Our 
prescription passes this test, since the one-point function we obtain is indeed independent 
of/. 
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4 Examples 



In this section we apply the general prescription to several concrete cases. The examples 
below are meant to illustrate the applicability of the real-time gauge/gravity prescription for 
computing time-ordered, retarded or Wightman correlation functions in a variety of back- 
grounds directly from the bulk theory. Notice that such correlation functions sometimes 
differ only by the form of their ie insertions (or other analyticity properties). Although the 
ie insertions are often set by hand, in the QFT they can be obtained from a formal deriva- 
tion which is briefly discussed in appendix lAl Any first-principles real-time gauge/gravity 
prescription should therefore also be able to correctly determine these ie insertions via bulk 
computations. The examples below show that our prescription indeed produces ie insertions 
that are always in agreement with field theory expectations (as described in appendix |A|) . 
which provides a nontrivial check of the prescription. 

4.1 Examples involving global AdSs 

For the examples in this subsection, we will consider a two dimensional CFT with a holo- 
graphic dual defined on a cylinder with metric 

ds^ = -dt^ + dcj)'^ (4.1.1) 

and a contour for the CFT which consists of the (j) circle times a path C in the complex 
time plane, with C being piecewise horizontal or vertical. The discussion can be extended 
straightforwardly to a CFT in d dimensions, but we restrict ourselves to d = 2 for now. 
Various possibilities for C are indicated on the left of Fig. [H We will compute the two-point 
function for operators inserted on the last two contours drawn in Fig. the first contour, 
with the indicated operator insertions, was discussed in [7]. 

As we mentioned in the section [21 the general idea is to 'fill in' the entire field theory 
contour with bulk spaces. In the case when all sources vanish along C, one can fill each 
horizontal segment of C with a segment of empty Lorentzian AdSs and each vertical segment 
with a segment of Euclidean AdSa. The metric on the Lorentzian segments is of the form: 

ds^ = -(r2 + l)de + 4^ + r^d4>^ (4.1.2) 

-|- 1 

and the Euclidean metric can be obtained by the replacement t = —ir. In this metric, 
surfaces of constant t or r have vanishing extrinsic curvature and the induced metric is 
independent of the signature of the spacetime metric. Therefore, the matching conditions 
for gravity are satisfied if we glue the Euclidean and Lorentzian segments together along such 
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Figure 4: On the left, various contours in the complex time plane. The vertical segments in 
the first two contours should be thought of as extending to infinity, yielding a vacuum state 
on the corner. The circles in the third contour should be identified; it is then a thermal 
contour. The crosses represent an example of the operator insertions we consider. On 
the right, we sketch the spacetimes consisting of piecewise Euclidean and Lorentzian AdSs 
that fill the given contours. One should impose matching conditions on the hypersurfaces 
between the segments. 

surfaces. The complete bulk solution Mc consisting of Lorentzian and Euclidean segments 
glued together along these constant t or r surfaces therefore satisfies all the conditions 
stated above and can be taken as a filling for the given contour. We have drawn such 
fillings schematically on the right of Fig. HI Note that these 'piecewise AdS' spacetimes 
may not be the only bulk solutions for the given class of contours, a point which we will 
come back to when we discuss black holes. 

By switching on boundary sources, we can perturb such backgrounds, with the provision 
that the matching conditions are satisfied for the perturbations as well. In the two examples 
below, we will add a massive scalar field in the bulk and compute a contour-time ordered 
two-point function of the dual operator. We will work in the approximation in which the 
scalar field is free and propagates without backreaction. 

4.1.1 Generalities 

Before considering specific contours, we first discuss some generalities regarding the solutions 
to the Klein-Gordon equation that are valid for each Lorentzian and Euclidean segment 
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Figure 5: The dots represent poles in complex frequency space and the curved line a Feyn- 
man contour that avoids them. 

separately. 

We start from the action 

S = ^j (f^^x - m^^>2) , (4.1.3) 

with d = 2 and the metric G^i/ given by (j4.1.2p . As usual, we have m? = A(A — 2) with 
A — 1 = / G {0, 1, 2, . . .}. As we already discussed in [7], the regular mode solutions to the 
Klein-Gordon equation are of the form 

e-''^'+'''^f{u;,\klr), (4.1.4) 

with 

f{u, k, r) = C^klil + r'^Y'^r^F{{uJ + k + l + l)/2, {uo + k + l- 0/2; A; + 1; -r^) , (4.1.5) 

where F is a hypergeometric function and C^ki is a normalization factor chosen such that 
the coefficient of the leading term equals 1. For large r, this solution behaves as 

/(w, fe, r) = r'-^ + . . . + r~'~^a(u;, k, l)[\n{r^) + /3(a;, k,l)] + ... (4.1.6) 
((c^ + A: + 1 - /)/2)f((cj - + 1 - l)/2)i 

-ij{{uj + A: + 1 + 0/2) - ^((-w + A: + 1 - l)/2) , (4.1.7) 

where {a)n = r(a + n)/T{a) is the Pochhammer symbol and iIj{x) = dlnT{x)/dx is the 
digamma function. In the expansion (j4.1.6p we omitted terms of lower powers of r and 
some terms polynomial in u and k (which would lead to contact terms in the 2-point 
function). 

A bulk-boundary propagator can be obtained by an integral over co and a sum over k of 
these modes. However, if the frequency equals 

uj = u^i^ = ±{2n + k + l + l), n G {0,1,2,...}, (4.1.8) 



with 



a{uj, k, 

(3(L0,k,l) 
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then the term a{uj, k, l)P{uj, k, I) in the radial expansion of the modes has a pole and therefore 
the modes become singular. To obtain a well-defined bulk-boundary propagator, one needs 
to specify a contour in w-space around these poles, for example the contour sketched in 
Fig. [5l Furthermore, the residues of these poles form exactly the normalizable modes, since 
the poles only occur at normalizable order in the radial expansion of the modes (indeed, 
the first terms in this expansion are always local so they cannot contain poles in frequency 
space). Let us denote them by g{u;nk, k,r): 



g{uJnk,k,r) = d) dujf{cvnk,k,r 

-l-l 



r 



a{LOnk,k,l)(^f dL0(3{u;,k,l)^ +... (4.1.9) 



= r~ ~ 4:TTia{uJnk, |A;|,0 + • ■ ■ , 
where the contour is defined as counterclockwise for w"^ and clockwise for w^^, so that 
ff(^nfe'I^U) =5(^^~fc,I^U) and q(u;+^„ l/cl,/) = a(a;~^, |A;|, /). 

These normalizable modes can be added at will to any solution without changing the 
asymptotics for large r. The most general solution (without specifying any initial or final 
data) therefore involves an arbitrary sum over these modes and is thus of the form 

Ht,cl>,r) = ^Yl 1^1 [ #e-^"^*-*>^'^'^-'^V(o)(i,0)/(c^,|A;|,r) 

(4.1.10) 

+ E E E -*+^^-<^ff(a;.., \k\ , r) , 

with so far arbitrary coefficients c^^ (provided the sum converges). For convenience, let us 
fix the contour C to be of the Feynman form sketched in Fig. [5l Any different contour can 
then be implemented by changing the normalizable modes. As we show in detail below, 
the initial and final data, so the other segments and matching conditions, will eventually 
completely fix the c^^. 

Below, we will often make use of the following observation. To the past of all the sources, 
the contour of the w-integral can be closed in the upper half of the complex frequency 
plane. The choice for a Feynman contour implies that we pick up the poles at the negative 
frequencies only, which we repeat are just the normalizable modes. The solution can then 
be fully written as a sum over normalizable modes, 

^ oo 

n=ofcez ^ (4.1.11) 

+ E E E de-^^-*+^'^5(^n., \klr) , 

± fceZn=0 
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Figure 6: The in- in contour we use to compute a Wightman function. We choose time 
coordinates that increase in the direction of the arrows. 

which is to be expected by completeness of the normahzable modes. Similarly, to the future 
of all the sources, we can deform the contour in the lower half plane and pick up the residues 
at the positive frequencies. 

Next, consider the solution on the Euclidean segments. One can obtain the mode so- 
lutions by a replacement of the form t = —it. We will set all sources to zero along the 
Euclidean segments, so the solutions there will always consist of normalizable modes only, 

oo 

$i,(r,,/.,r) = j;^^d±e-nV+^'=<^5(^„fc,|A:|,r), (4.1.12) 

± fcgZn=0 

with to be determined coefficients d^^. Note that, if a contour extends all the way to 
r oo, then we also require finiteness of the solution in this limit. This corresponds to the 
absence of any sources at this point. Such a condition directly implies that all the are 
zero, whereas the are still unconstrained. The converse statement holds for a contour 
extending to r — > — oo. 

This finishes the introduction of the solutions on the various segments; we can now 
consider specific contours and see how the matching conditions specify the coefficients of 
the normalizable modes for us. 

4.1.2 Wightman functions 

Our first example is the computation of a vacuum-to-vacuum two-point function using an 
in-in formalism. As explained in appendix [21 the in-in formalism in particular allows for 
the computation of Wightman functions directly from a path integral. In our case, we can 
do the same holographically. 

Let us therefore consider the contour sketched in Fig. given again in Fig. [6l It runs 
from ioo to 0, then to T (with T real and positive), then back to the origin and then to 
—ioo. As we outlined above and sketched on the right of Fig. [Jja, for such a contour we 
consider a filling that consists of two Lorentzian AdSs spacetimes between two Euclidean 
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AdSs caps. These four space(time)s will be denoted as Mj, with i running from to 3. 
We will use a subscript i also on other quantities to distinguish on which of the Mi they 
are defined. Sometimes, in order to avoid confusion with other subscripts, we will put this 
subscript in square brackets, writing for example C[j]. 

We can again split the contour-integrated action into the following combination: 

/O f-T r-2T roo 

dToLE{^o)+i dhLL{^i)-i dt2LL{^2)- dTsLEi^s), (4.1.13) 
-oo Jo JT Jo 

with the Lagrangians 

Li($) = i /" d'^xy^li-d^^d'''^ - m2$2) ^ 

(4.1.14) 

Le{<^) = d'^xy^\{d^<^>d''<S> + m2$2) . 

We use a real contour time coordinate on every segment Mj whose direction is indicated in 
Fig. [6l We glue the surface given by rs = to that given by t2 = 2T, and similarly the 
surfaces ti = T to t2 = T and ti = to tq = 0. 

A full list of matching conditions is now given by continuity of the fields, plus continuity 
of their derivatives with appropriate signs. These signs are easily found by equating the 
conjugate momenta obtained from functional differentiation of the on-shell actions. One 
obtains: 

-droM^o = 0) - idt.^iih = 0) = 
+idt,Mt2 = T) + idt,Mt2 = T) = (4.1.15) 
-idt,Mr2 = 2T) + dr^M^s = 0) = . 

Consider the case of a nonzero source '^(o)[i] only on the conformal boundary of Mi, so 
is given by (cf. (|4.1.10p ): 

$i(ti,<A,r) = [ duj [ /4e-*^(*^-*)+*'=(*-'^V(0)[i](i,<^)/('^,|A:|,r) 

oo 

+ E E E <ti]nk^~'^-''''"'^9{u^nk, \k\,r) , (4.1.16) 

± k&n=0 

with to be determined coefficients cj^j^;.- As in [7], we take the source to vanish near ti = Q 
and ti = T. By performing the w-integral, we write the solution as a sum over normalizable 
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modes in the vicinity of these hypersurfaces: 

^ oo 

n=OfcGZ 

oo 

+ E E E ^J]nfce"'""^*^^''^5K,fc, I r) , (4.1.17) 

± keZn=0 

^ oo 



n=0 fcGZ 

oo 



+ EEE4n.^"^^-*^^^''5('-n.,|A:|,r). 

± k&n=0 

Since there is no source on the other segments, the solutions $0) ^^2 and $3 are just sums 
over normahzable modes. For <J>2 we obtain: 

00 

$2(t2, = E E E 4nfce-*""^*^+*''^5Kfc, 1^1, r) . (4.1.18) 

± fcgZ n=0 

For Mq we can only allow for modes of the form e"''l'^l'^o, since tq extends to —00. Similarly, 
since ra — > 00 on M3, the modes there are of the form e"''^''^''. We thus find that 



^o(to, M=Y.Y. C[0]nfce'^"^""+*''^5(^nfc, \Kr 
'&3(r3, r) = J] J] C[3]„fce-'^n^^3+*'=<^5(^nfc, I 



00 

^ ^ ,Tj • 

keZn=0 



The matching conditions will now determine the c|^]„^ for us. Since the different modes 
g((u;„/fc, |/c|, r) are orthogonal (up to the symmetry g{u}^^,k,r) = g{uj~f^,k,r)), we can do 
the matching 'mode-wise', i.e. we can compare the coefficients of the various modes. For 
example, the first matching between Mi and Mq, which is ^i{ti = 0, (/>, r) = ^{){tq = 0, (/>, r), 
yields 

RO]nfc = ^'A(0)[1] fc, ^) + C[^]„fc + Cfi]„fc , (4.1.20) 

and the second matching condition, which is the first equation in ()4.1.15p . becomes 

- '^nfcClOjnfc - 4^<fc'/'(0)[l] k) - ^^+fcC+ - UJ~k(^[l]nk = • (4.1.21) 

Recalling that = —^jJ^k combining the two matching conditions, we find that 

c+„, = 0, (4.1.22) 
which is the statement that there are no positive frequencies to the past of the sources. 
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Similarly, from the matching conditions between M2 and M3, one deduces 

c+„, = 0, (4.1.23) 

so on M2 we can only allow for negative frequencies with respect to ^2- Then, from the 
matching condition between Mi and M2, we see that frequencies should be inverted on M2: 
positive frequencies on Mi become negative frequencies on M2 (with respect to t2) and vice 
versa. Therefore, on Mi there can only be positive frequencies close to ti = T. Indeed, 
working out the details results in 

which, combined with (j4.1.22p . completely fixes the c^]„^, on Mi to be zero. We thus 
obtain the usual Feynman prescription for the bulk-boundary propagator on Mi, which is 
reassuring: using the in-in instead of the in-out formalism should not have changed our 
result of [7] and indeed we found that it did not. 

The solution is now completely fixed and one may compute all of the c|^]„^ and the 
C[o]nfc ^iid C[3]„fc using the matching conditions. For the Wightman function, we will only 
be interested on the solution on M2 for a source on Mi, so we will only need the c^i^^f^- 
Equation (I4.1.23P already fixed half of them, and the first matching condition between Mi 
and M2 yields 

''[2\nk = ^Mmi^nk^ k)e-'<^^ . (4.1.25) 

With the solutions determined, consider the one-point functions. As we mentioned in 
section [3l the gluing of different solutions does not affect the usual prescription that the 
renormalized one-point function in the presence of sources is given by the renormalized 
radial conjugate momentum. For the case under consideration, we thus obtain (up to 
contact terms) 

(0[i](x)) = _^_i_(-5o + iSi - iS2 - 53 + Set) = -2Z</.(20h(x) , (4.1.26) 

with 4'{2l)[i] the term of order /--'-i in the radial expansion of 
We are in particular interested in the Wightman function: 

(0(x)0(x')) = (rcO[2](x)0[i](x')) 
i 6 

= 7^<50(o)[i](xO^^[^l^"^^ (4.1.27) 



-2li 
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With the solution ^2 we found above, we obtain 

00 „ 

(4.1.28) 

Using t = 2T — t2 and taking the functional derivative, we get: 

{Oix)0{x')) = ^ j;e-:.(*-*')+^'=(<^-<^')a(u;+ ,|fc|,0. (4.1.29) 

n>0 fcez 

This expression satisfies some standard checks that are expected for a Wightman function, 
namely it vanishes for a; < and the coefficients are real and positive definite, see appendix 
lAl Evaluating the summations, we find 

{0{x)0{x')) = ^'Jf''^ , (4.1.30) 

[cos(e — le) — cos((pjJ'+^ 

which has poles when t — ie\s real, so it is analytic in the lower half plane, also as expected. 
4.1.3 Thermal AdS 

Let us now consider the thermal contour indicated in Fig. U];. We again take ti to run from 
to T on Ml, t2 to run from T to 2T on M2, and r to run from to /3 on M3. The novelty 
here is that we glue the part with r = /3 to the surface with ti = in order to obtain a 
thermal state. The action splits into: 

j-T f2T 

i dtiLL[^i]-i dt2LL[^2]- dTLE[^z]- (4.1.31) 
JO Jt Jo 

Again, we consider a source living only on Mi and solve the Klein-Gordon equation. The 

general expressions for <I>i and ^2 (without specification of initial and final data) are exactly 

the same as before and are given by the equations (I4.1.16P and (|4.1.18p . We can also use 

the expressions (|4.1.17p for $1 when t ~ or t ~ T. For <^3, we may allow for both positive 

and negative frequencies and the most general purely normahzable solution is: 

00 

Mr, <^,r) = Y,Y.Il 4ink^^^'^^""^9{u^nk, \k\ , r) . (4.1.32) 

± keZ n=0 

The full list of matching conditions is now 

Mh = T) = $2(^2 = T) dt^ti =T) = -dtMt2 = T) 

$2(t2 = 2T) = $3(t = 0) idtMt2 = 2T) = dr^r = 0) 

$i(ii = 0) = $3(t = (3) -idt^ti = 0) = dr^T = 13) , 
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which, after some algebraic manipulations, results in 



'[l]nk 



1 

4^ 



'/'(0)[l] i^nk^k) 



1 



(4.1.33) 



The nonzero Cry, , directly enter into the two-point function and we get 



(rO[i](x')0[i](x));3 = / ^'^e—(*-*')+*'=('^-'^')a(a;,|fc|, 0/3(^,1^1,0 





(4.1.34) 



with the subscript f3 indicating the temperature. As in the 'free-field' approximation, we 
find the sum of the zero-temperature Feynman propagator and a heat-bath contribution. 
Also, notice the symmetry x x'. After rewriting the thermal contributions as geometric 
series, one readily finds that this expression becomes 



which satisfies the KMS condition, so it corresponds to a thermal two-point function indeed. 
It is a sum over images of the zero temperature result, reflecting the fact that Euclidean 
thermal AdS is obtained by identifications in the time direction of Euclidean global AdS. 

One can actually arrive more directly at (j4.1.35p by using the relation between thermal 
AdS and global AdS to first obtain the Euclidean correlator by a sum over images and 
then analytically continue to real-time. This was the way (j4.1.35p was obtained earlier in 
|42| . Of course, the ie insertions then have to be fixed by hand. The emphasis here is on 
the fact that we can unambiguously arrive at (j4.1.35p . including the correct ie insertions, 
by employing a Lorentzian signature gauge/gravity dictionary and without assuming any 
special properties of the background under consideration. 

4.2 Poincare coordinates 

For our next example we consider a CFT in d-dimensional Minkowski {Minkd) spacetime. 
As is well known, Minkowski spacetime is conformally isometric to an open region of the 
Einstein static universe, M x S'^~^ . Thus the correlators for the CFT in Minkowski space- 
time can be obtained from those of the Einstein universe (as we demonstrate for d = 2 in 
subsection I4.2.5P . Since boundary Weyl transformations are a specific class of bulk diffeo- 
morphisms |43^ I44j. a similar procedure can be done holographically. 

Nevertheless, it is still interesting to directly compute the correlators in Minkowski 
spacetime, not least because this is the typical background for most QFT computations. 




(4.1.35) 
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Furthermore, for a CFT on Mink^ to be exactly equivalent to the theory on M x S'^~^ 
the boundary conditions of all QFT fields at infinity of Mink^ must be the ones dictated 
by the theory on M x S'^~^. One may however wish to study the QFT on Mink^ with 
fields satisfying different boundary conditions at infinity. For example, the ground state 
of a conformally coupled scalar on M x S'^^^ has necessarily (0) = because of the 
curvature coupling of the scalar. The same theory on Minkd however allows for ground 
states with non- vanishing (0), since in this case the curvature coupling vanishes. In such 
cases the nonzero scalar vev spontaneously breaks conformal invariance. This is described 
in the bulk by domain wall spacetimes containing additional bulk fields capturing the vevs 
of gauge invariant operators. One can extend the methods described here to apply to the 
computation of real-time correlators along holographic RG flows, extending the Euclidean 
computations in |23^ I24j. but we shall not discuss this in detail here. 

Instead we will compute vacuum-to-vacuum amplitudes for the CFT without vevs. To 
this end, we consider the field theory path of Fig. [1^, but with M'^"^ as the spacelike part 
of the boundary manifold. We can compactify the entire contour by adding a single point, 
resulting in the boundary geometry shown in Fig. [7^. The Lorentzian segment is cut off 
at finite initial and final times t = ibT. Below we holographically compute a time-ordered 
two-point function for a CFT in this background. 

4.2.1 Bulk spacetime 

As before, the first step is to find a suitable bulk manifold that fills in the contour. 

We begin with the Lorentzian segment of the contour. In the absence of any sources 
and vevs, it is filled in with a segment of empty AdS^+i in Poincare coordinates: 



The Poincare coordinate system covers only a part of all of AdS^+i, as indicated in Fig. [TJa. 
We will however cut off the bulk manifold along the hypersurfaces t = itT and therefore we 
will not need the rest of the AdS^+i spacetime anyway. The Lorentzian segment with the 
above metric and — T < t < T will be referred to as Mi below. 

The two Euclidean segments can be filled with Euclidean AdSj^+i, whose metric can be 
obtained from ()4.2.ip by the replacement t = —it. We again need only a part of these spaces 
and cut off the Euclidean solutions along hypersurfaces of constant r. More precisely, we 
call Mq the Euclidean manifold with the metric 



dz^ - dt^ + dx^ 



(4.2.1) 



Z' 



.2 



ds^ 



dz"^ + dr^ + dx^ 



(4.2.2) 
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(a) 



(b) 



Figure 7: (a) The geometry used for the computation of the two-point function in Minkd- 
The Lorentzian manifold is cut off at shoes given by t = itT, to which the darker shaded 
Euchdean caps are glued, (b) The Poincare coordinate system covers only a part of AdS^+i- 
Both the global AdS time and the Poincare time run upward. The planes t = ±00 bound 
the coordinate system. 

and To < 0. Similarly, we take M2 the Euclidean manifold with the r2 > and the same 
metric (|4.2.2|) with the replacement tq ^ T2. 

Next, we glue the three components together by gluing the surface given by tq = on 
Mq to the surface t = —T on Mi, and the surface r2 = to the surface t = T on Mi. One 
may easily verify that the matching conditions for gravity are satisfied, since the induced 
metric on surfaces of constant t or r is the same and these surfaces are totally geodesic. We 
conclude that the combination of Mq, Mi and M2 satisfies all the holographic boundary 
data as well as all the matching conditions, and so it can serve as the background around 
which we study perturbations below. 

4.2.2 Solutions 

We will again obtain a time-ordered two-point function of a scalar operator of conformal 
weight A = I + with / G {0,1,2,...}, which is dual to a bulk scalar field of mass 
m? = A(A — d). As we did in the previous examples, we take the scalar field to propagate 
freely and without backreaction. 

On Ml the action for the scalar field is again (j4.1.3p . this time with the metric (j4.2.ip . 
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Solutions to the equations of motion satisfy the Klein-Gordon equation: 

z'^+^d^{z-'^+'^d^<^>) + z^Do^ _ = . (4.2.3) 
After separation of variables we find modes labeled by (w, k): 

For spacelike momenta = —uP' + A;^ > 0, these modes are unambiguously defined. For 
timelike momenta < 0, we have to consider possible branch cuts. First of all, we put the 
square root in defining q = just above the negative real axis. We indicate this by using 



-uj^ + k^-ie. (4.2.5) 

Second, Ki has a branch cut along the negative real axis, which is however unimportant 
since | arg((7ez)| < tt/2. Finally, Ii has no branch cut since I is an integer. 
To select the right solution on Mi, we should look at the asymptotics: 

^d/2-l 

z'^/'K,iqz^O)=m^ + ... 



I ^d/2+l 

r(/TT) W 



z^/^Kiiqz ^ oo) = ^ ^^e-"^ + . . . 



(4.2.6) 



z^/^Il{qz ^ oo) = J^le'^' + g-^^-^'+s)-'] + . . . 
y 27r(7 

For spacelike momenta, finiteness as z — > oo selects z'^^ Ki{qz) as the only correct solution. 
On the other hand, for timelike momenta no linear combination of the solutions remains 
finite as z — > oo, which means that any solution that does remain finite as z ^ cxd should 
be obtained as an infinite sum over the modes. Furthermore, from the asymptotics as 
z ^ 0, we find that the modes z'^/'^ Ki{q^z) ~ z'^l'^~^ correspond to sources on the conformal 
boundary, whereas the z'^Pli{q^z) ~ ^'^1'^+^ are the normalizable modes. 

For timelike momenta g^z = —i\q\z and we will henceforth rewrite the modified Bessel 
function of the first kind using lliy-i\q\z) = e-''^^/^Ji{\q\z). Although we could also have 
rewritten Ki{—i\q\z) = {iire^^'^/'^ /2)Hn\\q\z), we do not do so below, since z'^/'^ Ki{q^z) is 
needed for both spacelike and timelike momenta. We emphasize that Ki{q^z) is unambigu- 
ously defined for all real q^. 

Next, consider the manifolds Mq and M2, both with the Euclidean metric (j4.2.2p and 
< r < 00 and — 00 < r < 0, respectively. Although we will mainly work in position 
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space below, we will for completeness present the mode solutions here as well. First of all, 
the mode solutions on Mq and M2 are obtained by the usual substitution t — > —it in the 
Lorentzian modes (|4.2.4|) . Since we will not switch on any sources on these segments, the 
solutions on Mq and M2 need to be purely normalizable. As we just showed, this implies 
that only the modes z'^^'^Ji{\q\z) with < are allowed. Furthermore, since no operators 
are inserted at the points r ±00, we will also request finiteness of the solution in this 
limit. This implies a restriction to negative frequencies on Mq and to positive frequencies 
on M2. More explicitly, the solutions on these segments are built up from the modes 

el'^l"»+^^-V/V;(|g|z) on Mo, 

(4.2.7) 

e-l'^l"2+^'=-V/2j/(|g|z) onM2, 

with — oj^ + k^<0. Since the individual modes diverge as 2; — > cxd, we should again sum an 
infinite number of these modes in order to get a solution that vanishes also at this point. 

4.2.3 Bulk-boundary propagator 

The next step is to compute a bulk-boundary propagator, which we denote by X{t,x,z). 
We will consider the propagator for a source on the conformal boundary of Mi only. Let us 
first investigate the solution on Mi. Inspired by the Euclidean bulk-boundary propagator, 
we may try: 

X,{t,x,z) = ^l^d^l dke-^^'+''-^^^!^z^/'Ki{q,z). (4.2.8) 

The ie-prescription is equivalent to a Feynman contour C in the w-plane around the branch 
cuts which we show in Fig. [8l The expression (j4.2.8p is not obviously convergent as z ^ 00. 
However, we can perform the Fourier transform by closing and deforming the contour. The 
ie-prescription tells us which branch cuts we pick up and the corresponding position-space 
expression is equal to 

Xi(t,x,z) = ir(/)r(/ + -)^-i ^2^^^ (4-2-9) 

which clearly converges for large z. 

As in the previous section, this bulk-boundary propagator is not unique without im- 
posing initial and final conditions. Indeed, one may always add a normalizable solution, 
which we will denote as Y{t, x, z). Notice that we know that normalizable solutions on Mi 
exist from our discussion of the previous section, where we used global coordinates. In the 
previous subsection we found that Y{t, x, z) must be a linear combination of the modes 
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Figure 8: The contour around the branch cuts (wavy hues) in the complex frequency plane 
used to define a bulk-boundary propagator. 

z'^/'^Il{q^z) with < 0, which we write as 

Yiit,x,z) = jdu; I dke~''^'+''^-^e{-q^)c^^^{LJ,k)z''I^Ji{\q\z) , (4.2.10) 

with further constrains on C[i](a;,A;) by requesting finiteness for z ^ oo that we will not 
work out here. To reiterate, without initial or final conditions such normalizable solutions 
can be added at will to our suggested bulk-boundary propagator ()4.2.8p . so the normalizable 
solutions parametrize the ambiguity in the bulk-boundary propagator. In particular, any 
different ze-prescription than the one we fixed above can be implemented by changing these 
C[i] {i^,k). 

4.2.4 Matching 

With the solutions on Mi specified, let us now discuss the matching. We will show that the 
matching conditions imply that Xi(t, z) is the right bulk-boundary propagator and that 
no normalizable solution can be added since Yi(t,x, z) can never be matched to a regular 
and normalizable solution on Mq and M2. 

We begin with the matching conditions between Mq, Mi and M2: 

$i(ti = T,x,z) =$2(r = 0,x,z) idt^^i{ti = T, f, z) + a^^>2(T = 0, x, z) =0 

^i{ti = -T,x,z) = $2(r = 0,x,z) -idt^^iih = T,x,z) - dr^2{T = 0,x,z)=0. 

(4.2.11) 

Let us now show that we can find solutions Xq and X2 on Mq and M2 that can be matched 
to Xi. This is straightforward in position space, where we can verify that the position-space 
expressions 

Xo{To,X,z)=iT{l)r{l + -)TT~^ 



2' {-{-T - iror + ^ + z^ + ieY^''/^ ' ^ 12) 

X2iT2,X,z) =iT{l)T{l + ^)TT-i—— . x2^'-^^ 2^ ■ \l+d/2 ' 
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satisfy the equations of motion on all of Mq,M2 and are normalizable. Furthermore, they 
actually satisfy the matching conditions as well. To see this, notice that the (+ie)-insertions 
in the denominators of (|4.2.12|) are not necessary for nonzero r, but they are necessary to 
ensure that ()4.2.12p are well-defined (distributions) on the initial and final hypersurfaces 
given by To = and T2 = 0. With the given ie-insertions, we can compare (14.2. 12p to (I4.2.9P 
and one readily verifies that the matching conditions are satisfied. 

Let us now show that one could not have picked any other ze- insertions {—ie, +iet, 
etc.) on the Lorentzian side. If we would have done so, the matching conditions would 
directly dictate a corresponding change in (j4.2.12p . However, such a change in the Euclidean 
solutions is not allowed, because any other ze-insertion in (j4.2.12p would give a singularity 
in either Xq or in X2. For example, if we would replace the +i€ with —ie on M2, then X2 
would be singular at T2 = e/2T, around the point given by + = and thus this 
solution should be discarded. We conclude that the ze-insertion in (j4.2.9p is the only one 
that moves the singularity everywhere away from the contour. 

It remains to show that (j4.2.9p is the indeed the unique bulk-boundary propagator 
by demonstrating that there are no matching Euclidean solutions for the normalizable 
solution (j4.2.10p . Using the normalizable modes we found above, the solution on Mq should 
necessarily be of the form 

Yo{to,x,z) = j du j (ifeel"l^«+*^--"^0(-g2)^^^^(|^|^^)^d/2_^^(|^|^)^ ^42.13) 

for some coefficients C[o](|u;|, /c). A similar expression holds for the solution on M2: 

Y2{t2,x,z) = j du j dfce-H-2+*fc-0(_^2)^^^j(|^|^^)^d/2_^^(|^|^)^ ^4_2.14) 

Consider now the matching conditions, for example the continuity condition between Mi 
and M2: 

Yi{T,x,z) = Y2{to = ^). (4.2.15) 
Although this is an equality between two integrals, the modes z'^^'^ Ji{\q\z) are orthogonal, 

dzz-^Jl{\q\z)Ji{\q'\z) = cd{\q\ - \q'\) , (4.2.16) 

with c a constant. We can therefore equate the integrands (up to w <-> —lo), which results 
in 

C[i]{uJ,k) + C[i]{-uj,k) = C[o]{\uj\,k) . (4.2.17) 

The other matching conditions can be imposed in a similar way and they ultimately deter- 
mine C[i](w, k) = C[o](u;, k) = C[2](<^, k) = 0. There is thus no normalizable solution and the 
bulk-boundary propagator X(t,x, z) is unique. 
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4.2.5 Two-point function 

As for the computation of the time-ordered two-point function, the only difference with 
the Euchdean case are the ie-insertions in the bulk-boundary propagator, which in Fourier 
space corresponds to the replacement q — > g^. Just as for AdS^f+i in global coordinates, 
these enter directly in the two-point function which, up to contact terms, is then given by: 

{TO{q)0{-q)) = ^Z^qflogq,. (4.2.18) 

The ze-insertion corresponds again to the Feynman contour of Fig. [H] around the branch 
cuts, signifying time-ordering indeed. In position space, we find 

2ir{l + d/2) 1 



d ' 

2 



(4.2.19) 



and the ie-insertion agrees with 



Normalization 

Let us compare the normalization of the time-ordered two-point function on the cylinder 
with that on two-dimensional Minkowski space. We start from the time-ordered two-point 
function on the cylinder: 



2 I Jj,2 



[cos(i — iet) — cos{(j))y 
and apply the coordinate transformation t = u — v, (j) = u + v, after which we find 

{TO{x)Om = —- /\ . ,.+1 , ds' = Adudv , 

[sin(n — 17]) sm[v + zryjj'^^ 

where now rj = e{u — v). We then Weyl transform and use covariance of the two-point 
function: 

,^.^/ N.^/ XV 2/2 /vr o dudv 

{TO{x)Om = J—, . , / , , . , ds^ - 



[tan(n — zry) tan(7; -|- ir/)]'+^ ' cos2(m) cos2(u) ' 

where we should remember that the two-point function is multiplied by two Weyl factors; 
one evaluated at x and one at 0. For — 7r/2 < u,v < 7r/2, we can rewrite the denominator 
using 

tan(n — irj) tan(f -|- irj) = tan(ti) tan(u) -|- ie{u — v)[tan{u) — tan{v)] = tan(u) tan(u) -|- ie' 
with e' positive and constant. Finally, using x + y = tan(ti) and x — y = tan(t;), we obtain 
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and the normalization is indeed the same as in (|4.2.19|) evaluated at d = 2. 



4.3 Higher-point correlation functions 

In this subsection, we briefly discuss how real-time higher-point correlation functions can 
be computed with our prescription. We take an interacting scalar field with potential 

2 3 
so that the equation of motion becomes: 

- m'^<^> - = . (4.3.2) 

This equation can be solved perturbatively. We first compute the sequence: 



{0} 




= 0, 


{1} 


- m2$^i} 


= A'^^o} 


'{2} 


- m^^{2} 


= ^^{1} 



(4.3.3) 



where ${o} satisfies the radial boundary data and <I>|j} with i > 1 vanish asymptotically. 
The full solution is then obtained as: 

$ = ${0} + ${1} + ^>{2} + • • ■ (4.3.4) 

To compute the series we need to compute the bulk-bulk propagator Z. This propa- 
gator satisfies 

aGZ{x,x') = ^<5''+' {x - x') (4.3.5) 
and vanishes asymptotically. In terms of Z(x,x'), we find: 

<I>{i+i}(x) = / d''+^x'^Z{x,x')^^}{x'). (4.3.6) 

J M 

In our case, the bulk manifold M splits into multiple parts and we need to integrate the bulk- 
bulk propagator against <I>|j} on the various segments. The bulk-bulk propagator therefore 
also splits in multiple components depending the segment that x and x' lie on. We will 
indicate this by a subscript. For example, Z[i2](x,x') denotes the bulk-bulk propagator 
with X on Ml and x' on M2. Equation ()4.3.6p then becomes: 

^[i]{i+i}{x) = Y.f d'^+ix'V^Z[,fc](x,x')$[fc]{.}(^0 , (4.3.7) 
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with the sum over all of the components Mj.. Of course, Z\^j^ is homogeneous on Mj if 
j 7^ k. Also, we wih exphcitly see below that x') = x). 

Let us now find this matrix of bulk-bulk propagators. These bulk-bulk propagators need 
to satisfy the matching conditions, since then so will all the ^{i] and consequently also <I>. 
(Our derivation of the matching conditions for a scalar field in section HTTl was independent 
of the potential V^[<I>], so the matching conditions are unchanged by the interaction terms.) 
For concreteness, consider the bulk spacetime of the previous subsection, with a Lorentzian 
segment Mi sandwiched between two Euclidean segments Mq and M2. The matching 
conditions become important when we move x from, say Mi to AIq while keeping x' fixed. 
For example, we get 

Z[ii\{ti = -T,x,z;t'i,x' ,z) = ZyQ^{TQ = Q,x,z;t'i,x' ,z') , 

(4.3.8) 

-i5t Z[ii](ti = -T,f, z;4,x',2;') - 5rZ[oi](ro = 0, x, z; 4, f', z') = 0, 
just as in (j4.2.1ip . and all the other matching conditions are similar. 

The uniqueness of the bulk-bulk propagator is clear from the previous section, where 
we showed that there is no normalizable homogeneous solution that satisfies the matching 
conditions. As for existence, the bulk-bulk propagator for Lorentzian AdS in Poincare 
coordinates is already known, see for example [5] and the references therein, where one may 
find that 

Z^ii^{x,x') = Z[iii], (4.3.9) 
with ^11 an AdS-invariant function. 



^ {z-z'f -(t-t'Y + {x-x'Y + ie 
= ,„„,r'l^L fl - A - I + l; [1 - fnl^) . (4.3.11) 



and Z given by 



vrd/2r(A - §)r(2A - d) V "^V ^2' 2 ' 2 
with -F(o, h; c; z) the hypergeometric function. This solution is regular except when — > 0. 
Analytically continuing t, t' to Mq or M2 by the replacement t = —T—itq or t = T—iT2 yields 
other ^ij, and the ie- insertions again ensure that these satisfy the matching conditions 
when either x or x' moves from one segment to the other. So if we define 

%](x,x') = Zfe,-], (4.3.12) 

then the various Zjjj] satisfy (j4.3.5p . the matching conditions, and vanish asymptotically. 
Therefore, the full matrix of bulk-bulk propagators can be obtained by this analytic con- 
tinuation. Just as for the bulk-boundary propagator, the matching conditions uniquely fix 
the ie- insertions to be those in equation (j4.3.10p . 
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Figure 9: The contour we use for the black hole. The circles should be identified. 

Again, these ie-insertions enter directly into the higher-point correlation functions. 
These are obtained as usual by further functional differentiation of the renormalized one- 
point function. For example, for a time-ordered vacuum-to-vacuum three-point function 
with all three arguments on Mi we obtain 



bulk solution, the procedure to obtain these correlation functions is therefore just as for 
Euclidean metrics, except for the replacement t — > t — zet (so — — > — -|- ie). 

4.4 Stationary black holes 

The thermal contour drawn in Fig. [Tt admits another possible bulk solution, which cor- 
responds to an eternal black hole. In this section, we will use this filling to compute the 
time-ordered two-point function for an operator dual to a free scalar field moving in the 
black hole background. We will again work in = 2, so the bulk spacetime is the static 
three-dimensional BTZ black hole. The rotating black hole will be discussed in the next 
subsection. 

Below, we will actually use the deformed contour of Fig. [9] rather than the contour of 
Fig. [It- As we will see shortly, this has the advantage of 'opening up' the second boundary 
of the black hole spacetime as well. In the next subsection, we describe a bulk manifold 
that fills in this deformed contour. Afterwards, we proceed by switching on a scalar field 
and holographically compute correlation functions. 

4.4.1 Bulk spacetime 

Consider the eternal Lorentzian massive non-rotating BTZ black hole, whose Penrose dia- 
gram is given in Fig. [TOk . The black hole splits into four parts, which we denote by L, R, 




(4.3.13) 



with (p(^2A-d) the coefficient of the normalizable mode (of order z ) in the z-expansion of 
and the source (/>(o) should be set to zero after the functional differentiation. Given the 
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(a) (b) 

Figure 10: (a) The Penrose diagram for the eternal BTZ black hole. The arrows indicate 
the direction of time. In the diagram, every point represents a circle. The horizons, which 
are the solid diagonal lines, separate the spacetime in four regions labelled by L, R, F and 
P. (b) We cut off the spacetime along the dotted lines and keep the part in between them. 



F and P. On either part the metric is 

ds^ = -{r^ -rl)dt^ + —f^-^ +r^d(l)'\ (4.4.1) 



If necessary, we will use a subscript like L or R to indicate the corresponding part of the 
spacetime. Notice that time runs backward on R. The mass and temperature of the black 
hole are given by 

2 

M = ^ , T = ^ . (4.4.2) 

8G3 ' 27r ^ ' 

(Recall that we set the AdS radius to one, £^ = 1.) To simplify the notation, we make the 
coordinate transformation 

t = — , r = r'r+ , </>=—, (4.4.3) 
after which the metric reads 

(r^ — 1) 

where we have dropped the primes. Note that the periodicity of (j) has now changed to 

(j)r^ + 27rr+ . (4.4.5) 

At the very end of the computation we will return to standard conventions. 

To use this spacetime as a filling for (a part of) the contour of Fig. [9l we first have to 
cut it off along an initial slice, which we take to be the ti = = Q slice, as well as a 
final slice, which we choose to be the rp = f slice, with r < 1 a constant. These segments 
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(a) 



(b) 



Figure 11: (a) The Euclidean BTZ black hole, where again every point represents a circle. 
EucHdean time r runs as indicated, (b) We cut off the spacetime along the dotted hne given 
by r = and r = vr and keep the lower part. 

are the (blue) dotted and the (red) dashed lines of Fig. [TOb . respectively. As is shown in 
Fig. [TOb . we keep the segment in between these surfaces. Notice that > but t/j, < on 
this segment. We will need two copies of the segment, which we denote by Mi and M2. 
Next, consider the Euclidean solution with the metric 

ds^ = (r^ - l)dT^ + + r^dcP^ (4.4.6) 
(r^ — 1) 

and with periodicities 

T~r + 27r, (/)~0 + 27rr+. (4.4.7) 

Topologically, this solution is D2 x S^, with D2 a two-dimensional disk and the is 
parametrized by (p. As shown in Fig. \TT\ we will cut it in half along the hypersurface given 
by r = and r = tt, and keep the part given by < r < vr. We will again need two copies 
of this part, which we denote as Mq and M3. In Fig. [TTb . we have drawn these spacetimes 
as half a disk. 

We now glue the four manifolds together as shown in Fig. [121 Notice that Mq is glued 
to Ml such that the part with r = is glued to the part with = 0, and the part with 
r = TT is glued to the part with = 0. The same holds for the gluing between M2 and M^. 

Let us verify that the matching conditions for gravity are satisfied. First of all, the fact 
that Ml and M2 are identical means that the matching conditions for gravity are trivially 
satisfied along their gluing surface, which is the (red) dashed line in Fig. [T2j In fact, we 
could have glued Mi and M2 along any spacelike bulk hypersurface extending all the way 
to the two radial boundaries (and disjoint from the surfaces tL = tji = 0), and the matching 
conditions would still be satisfied. 



56 



Figure 12: The four components Mq, Mi , M2 and M3 are glued together to create a manifold 
that fills the contour of Fig. [9l The direction of the various time coordinates is the same as 
in Fig. M 

For the matching between the Euclidean and the Lorentzian segments, one may directly 
see from the metrics (I4.4.4p and (I4.4.6P that any surface of constant t or r has the same 
induced metric. One may also use reflection and translation symmetry to find that the 
extrinsic curvature of such slices must vanish. Therefore, the matching conditions for gravity 
are satisfied for this gluing, too. Finally, by passing to a coordinate system that is regular 
everywhere at the gluing surface, one may verify that there are no problems at the coordinate 
singularity at r = 1, either. 

Let us now turn to the matching conditions for a scalar field. The overall action (j4.1.3p 
can be split into a separate piece for each segment: 

iSi - iS2 -So -S3. (4.4.8) 

Continuity and the saddle-point approximation for the combination of actions (j4.4.8p de- 
termines the matching conditions to be: 
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(4.4.9) 



Incidentally, one may have wondered why the second set of horizontal line segments in Fig. [9] 
points to the left rather than to the right. This can be seen from the matching conditions 



57 



(|4.4.9|) : as one may verify they correspond to continuity in the complex time plane only 
if the contour has the shape of Fig. [H One may also verify that a replacement —tR2 
has no effect on the shape of the contour. 

4.4.2 Mode solutions 

We can now turn to the computation of two-point functions. We start by finding mode 
solutions to the Klein-Gordon equation, 

□G$-m2$ = 0, (4.4.10) 

on the various components. As usual, m? = A(A — 2) and we assume A = 1 + ^ with 
I G {0, 1,2, . . .}. In Lorentzian signature, we find two possible solutions, which we denote 
by V'i, 

= e-^'^*+*'='^/(±w, k, r) , (4.4.11) 

with a radial part given by 

1 \ iLo/2 



f{uj, k, r) = C^kl 1 ^ r 



-l-i 



X F{^{lo -k) + i(l + /), '-{u: + k) + 1(1 + l);iu: + 1; 1 - 1) , (4.4.12) 

with -F(a, h] c; z) a hypergeometric function and 

^ r(i(^ + fc) + i(l + 0)r(i(^-A:) + l(l + /)) 

= r(z^ + i)r(0 ^^-^-^^^ 

chosen such that the coefficient of the leading behavior of /(iw. A;, r) as r ^ oo equals one. 
The asymptotic expansion of the modes is given by 

= ^-i^t+ik^ (^^i-i + . . . + a{±uj, k, l)r~^~^ [ln(r2) + /?(±u;, A:, /)] + ...) , (4.4.14) 

with 



, n _ ( ^ Ah{^ + k) + \{l-l)W,{u-k) + \{l-l)), 

^ ' '^"^ /!(/-l)! 

/3(t^, k, I) = -i^C-ico + k) + i(l + /)) - V(^(u; -k) + i(l + /)) + local . 



(4.4.15) 



where the local terms we omitted from P{uj, k, I) originate from the expansion of the prefactor 
(1 — up to the relevant order. Such local terms lead to contact terms in the two- 

point function and will be omitted. The similarity between the modes (|4.4.12p and (j4.1.5p 
is not accidental: one may verify that the backgrounds with the metrics (j4.4.4p and (j4.1.2p 
are related by analytic continuation in complex (t, r, (/))-space, and so are the corresponding 
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mode solutions in these backgrounds. Since the behavior of the modes in the interior of the 
spacetime is different, we will not use this fact here. 

Near the horizon both modes oscillate infinitely rapidly. To see this, we transform to 
Poincare coordinates, given by 



tanh(t) 



y_ 

X 



2<h 2 2,2 

e ^ = X — y + z 



which brings the metric to the form 



^{dx^ - dy"^ +dz'^). 



(4.4.16) 



(4.4.17) 



With this definition the future and past horizons on the L quadrant are mapped to x = —y 
and X = y, respectively. Taking the near-horizon limit x it y ^ 0, we find 



i}± = C^ki exp ( ^ y ln(a; ± y)^ + i(A; ^ a;)(/> ) (1 + . . .) . 



(4.4.18) 



We can create modes that are well-defined almost everywhere on the Lorentzian segments 
via analytic continuation across the horizons, in the way specified by Unruh |21l I46j . see 
also [13] . Depending on whether we analytically continue from L to R via the lower or the 
upper half of the complex y plane, an extra factor of e'^^ or e~'^'^ should be added to the 
L mode to produce an R mode. Since this is the case for both and ip-, we find four 
different combinations: 



-ILOt+i. 



on L 



f{iO, k, r) on R 



-iLut-{-i. 



on L 



-iujt+ik<j)—TTU! 



f{uj, k, r) on R 



(4.4.19) 



''"f'f{-uj,k,r) 



on L 



-iLOt+ik(j>+TTU) 



f{—uJ, k, r) on R 



-ILOt+i. 



''I>f{-u;,k,r) 



on L 



-iujt-\-ik4>—7ruj 



f{—LO, k, r) on R . 

These modes form a complete set both on L and on R, and can thus be used to decompose 
any solution. In particular, solutions that are regular at the horizons can be obtained as an 
infinite sum over these modes. 

Finally, on the Euclidean solutions Mq and M3, the mode solutions are as usual obtained 
by the replacement t — > —ir in the ■0±- In this case, there is no need for an analytic 
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continuation, and we find two rather than four solutions, which we denote by 

cj)^ = e^^+'^^f{±uj, k, r) . (4.4.20) 



Going through the same arguments as before, one finds that these modes also oscillate 
infinitely fast near the horizon. 

4.4.3 No normalizable solution 

Let us now show the absence of a normalizable solution satisfying the matching conditions. 
This would imply uniqueness of any solution satisfying given radial boundary data. 
We begin on Mi where we write 

Yi = ^ / (iw(c[i] ++(/>++ + C[i]+_(/)+_ + C[i]_+0_+ + C[i]__(/)__) , (4.4.21) 
k 

with the C[i]-|-± some functions of lo and k. Notice that the sum is over r+Zc G Z to comply 
with the periodicity (|4.4.5p . The solution looks different in the various regions. Using 
(|4.4.19p we obtain 

yi,L = /d^e— *+*^*[(c[i]++ + C[i]+_)/(a;) + (c[i]„+ + C[i]„„)/(-^)], 

k 

Yi,R = 5^y"d..e— *+^'=*[(c[i]++e-- + C[i]+_e---)/(a;) 

k 

+(c[i]_+e-- + C[i]„_e--)/(-u;)] , (4.4.22) 

where here and below we suppress the k,r arguments from f{uj,k,r) for notational sim- 
plicity. By substituting the asymptotic behavior (14.4. 14p of the modes, we find that Yi is 
normalizable on both L and R if 

C[i]++ + C[i]„„ + C[i]+„ + C[i]„+ = , ^^^^^^ 
(c[i]++ + c[i]_+)e-- + (c[i]+_ + C[i]__)e--- = . 
Similarly, on M2 we consider 

Y2 = ^ duj{ci2]++(p++ + C[2]+_(/)+- + C[2]„+<?!'_+ + C[2]„_(/)— -) (4.4.24) 
fe 

and the same argument as above leads to the the same conditions (|4.4.23p but with C[i]-|--i- 
replaced by C[2]±±- Besides satisfying the same radial boundary data, the matching con- 
ditions between Mi and M2 imply that Yi and Y2 also have the same initial data on the 
matching surface. Since the solution on either Mi or M2 is uniquely specified by boundary 
and initial data, we find that C[2]±± = C[i]-|--|-. 
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On the Euclidean parts Mq and M3 the solution should be a linear combination of the 
Euclidean modes ()4.4.20p . We write it as 

^ (4.4.25) 

k 

As for Yi and Y2, the demand for normalizability implies 

C[o]+ + C[o]_ = 0, C[3]+ + C[3]_ = . (4.4.26) 

We now impose the matching conditions (j4.4.9p between the Euclidean and the Lorentzian 
solution. Using the orthogonality of the normalizable modes, this leads to algebraic relations 
between the individual coefficients c\^i]± and c^jT^±±. In particular, the matching conditions 
between Mq and Mi determine 

C[i]+- = 0, (4.4.27) 

while those between M2 and M3 fix 

C[i]„+ = 0. (4.4.28) 
Using (I4.4.23P we conclude that all the C[i]-|-± = and thus no normalizable solution exists. 

4.4.4 Bulk-boundary propagator 

We will now find the bulk-boundary propagator for a delta-function source at (t, (j)) on the 
L part of Ml. Since we have just shown the absence of any normalizable solution, any bulk- 
boundary propagator that satisfies the matching conditions is guaranteed to be unique. Let 
us therefore make an educated guess and consider a solution Xi on Mi that contains only 
the modes 0++ and (f) : 

^1 = ^ / + , (4-4.29) 

with new coefficients a[i]±± which are to-be determined functions of uj and k. Again, to 
comply with the periodicity of (f) given in (|4.4.5p . we need r+A; G Z as well as the extra 
prefactor of 1/^+ to normalize the boundary delta function. Notice that we already split 

off a factor e*'^*-*'=<^ from and ajij . On the two regions R and L, our ansatz takes 

the following form: 

^i,L = ^E / 'i'^e--(*-*>^'=(^-<^)[a[i]++/(u;) + a[i]__/(-u;)], 

/ ' . , (4.4.30) 

Xi^R = / 'i'^e— (*-*)+^'=(*-'^)[a[i]++e--/(^) - a[i]_„e---/(-u;)] . 
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As we mentioned above, we put a delta-function source on {t, (j)) on the conformal boundary 
of L and no sources on the conformal boundary of R. Substituting the asymptotics (j4.4.14p , 
such boundary conditions for Xi lead to 

+ «[!]- = 1^ (4.4.31) 

Notice that these conditions already fix the solution on Mi to be: 

-1 e^'^'^ 

«[!]++ = ^2.u _ 1 = ^2nu. _ 1 • (4-4-32) 

In passing, we mention that it is not manifest that Xi is finite at the horizons. To check 
this, one substitutes the near-horizon expansion (j4.4.18p of the modes and then computes 
the w-integral by contour deformation. One finds that an ie-insertion is necessary to ensure 
convergence and to regulate the lightcone singularity. (A subtle point is that the 

and a[i] both have a pole at w = 0, but the residues cancel each other so the contour 

can be freely deformed around this singularity.) The sum over k can be computed using 
similar methods as we employ for the two-point function below and the computation then 
shows that after the ie insertion Xi is regular at all the horizons indeed. Notice that the 
light-cone singularity is expected; we found a similar singularity when we wrote down the 
position-space expression (j4.2.9p in Poincare coordinates. It can be removed by integrating 
the delta function on the boundary against a smooth source. 

Let us now verify that we can find normalizable solutions on Mq, M2 and M3 such that 
the matching conditions are satisfied, so that Xi is indeed the bulk-boundary propagator 
on Ml. We start with the matching solution Xq on Mq. It should be a linear combination 
of the modes 4)± , 

^0 = E / du;e-^^"''t>e'^'-"^haio]+f{^) + «[o]-/(-^)) • (4.4.33) 

Let us consider the following coefficients: 

a[o]+ = a[i]++ , a[o]_ = -0[i]++ = a[i]__ - 1 , (4.4.34) 

with a[i]±-|- as given above. As one may directly verify by substituting the asymptotic 
behavior (j4.4.14p (now with t = —it), the solution Xq is normalizable since a[Q]_|_-|-ajo]_ = 0. 
Notice furthermore that < r < vr on AIq. Therefore, despite the factor e^'^ , the w-integral 
is still convergent along the real axis on Mq, because ~ e~'^'^^ for large positive uj. 

To verify that the matching conditions are satisfied between Mq and Mi, notice that 
the difference between the Euclidean and the Lorentzian solution on L, 

Xi^L{t = 0) - Xo(t = 0) = Yl I dLoe''^^-'''^(p— = (4.4.35) 
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since t > 0, so one can deform the contour in the upper half of the complex w-plane where 

(j) has no poles even at normalizable order. (Actually, near the horizon, it is the oscillating 

behavior of radial part of the modes that determines where to deform the contour to. Since 
this is still the upper half plane, the difference vanishes there as well.) A similar argument 
shows that the second matching condition on L as well as the both matching conditions on 
R are also satisfied. 

Next we consider the solution on M2, 

^ k 

(4.4.36) 

Since the radial boundary data on Mi and M2 are now different, we cannot use the argument 
used earlier for the normalizable solution Y to argue that a[2]±± is the same as a[i]±± and 
we have to compute a\-2\±±. 

To begin with, notice that the matching between Mi and M2 takes places on the F 
component of the black hole as indicated in Fig. [TOj Starting from the L quadrant one must 
cross the future horizon but not the past horizon to arrive at the F quadrant. Therefore, the 

modes and 0-| , which become singular at the future horizon, acquire an additional 

factor of e^'^^ as we move from L to F. However, the modes 4>-± become singular only at 
the past horizon and do not get such a factor. These factors should be included both on 
Ml and M2 and show up in the matching conditions: 

a[2]+_e-- + a[2]++e-- = a[i]+_e-- + a[i]++e-- = -^^^ , 

(4.4.37) 

a[2]-+ + a[2]— = + «[!]— = ^2^^ _ ^ ■ 

These two equations, together with those arising from normalizability on both sides of M2, 
completely fix the a\^2]±± to be: 

a[2]++ = 0, «[2]+- = ^2^0. _ I ' «[2]-+ = 0, «[2]— = ^2-.u _ ^ ' (4-4.38) 

and we have found a normalizable solution X2 on M2 that matches to the solution Xi on 
Ml. 

Finally, we need to verify that we can obtain a normalizable solution X3 on M3 that 
matches to X2. Since X2, in contrast with Xi, is already fully normalizable, X^ can be 
easily obtained by a simple analytic continuation of the solution on X2- Just as for Mq, 
one may again verify that the w-integral in X3 is convergent along the real axis, that X-^ is 
normalizable and that the matching conditions are satisfied. 
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Thus, the bulk-boundary propagator Xi can be matched to normahzable solutions on all 
segments. Since there are no solutions that are everywhere normahzable, we have obtained 
the bulk-boundary propagator for the black hole filling of the contour of Fig. O The 
same bulk-boundary propagator was actually written down in p^], where it was obtained 
by imposing boundary conditions at the horizon which are natural from considerations of 
quantum field theory in curved space |46]. We have now derived that this is indeed the 
correct bulk-boundary propagator for the real-time gauge/gravity dictionary. 

4.4.5 Two-point functions 

We are mostly interested in the time-ordered and Wightman function for real times. By 
looking at Fig. [U we find that we need operator insertions on the L component of either 
Ml or M2, because these segments lie along the real time axis. To simplify the notation we 
omit the subscript L, which should be understood in all formulas in this subsection. 

The one-point function in the presence of sources is again just the normahzable compo- 
nent (/)(2i) of the bulk-boundary propagator, times a factor —21 which is fixed by holographic 
renormalization, see section [3j Completely analogous to the analysis in section UJl this nor- 
mahzable component (p(^2i) can be read off by substituting (|4.4.14p in the solution Xi^l or 
X2^L- The two-point function computation is again completely analogous, and we find 

[a[i]^_i_a(ti;, k, l)l3{io, k, I) + a^i] a{—uj, k, /)/3(— k, I)] . (4.4.39) 

We recognize the structure of a time ordered propagator at finite temperature |47] . Such a 
propagator is of the form 

A(a;, k) = -n(a;)AA(w, k) + {1 + n(a;))A/j(w, k) , (4.4.40) 

with n{io) the Bose-Einstein distribution, 

-(^) = ' (4-4.41) 

and Ar and are the retarded and advanced thermal propagators, which should be 
analytic functions in the respectively upper and the lower half of the complex uj plane, see 
appendix [Al Since /? = 27r in our coordinates, we find a[i]++ = —n{uj). The structure of 
(j4.4.39p thus agrees with expectations. 

To obtain a position-space expression, choose t' = and t > 0. This allows us to perform 
the (j-integral by deforming the contour to the lower half plane and picking up the poles. 
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These poles come from (3{—uj,k,l) and from the a+_(_ and the a The former have poles 

at the quasinormal frequencies, 

uj = 00^^ = -i{2n + 1 + 1) ±k, (4.4.42) 

and the latter have poles at u; = —im with m G {1, 2, . . .} (the apparent pole at w = in 

a++ and a has zero residue). 

Afterwards, we compute the sum over k as follows. We first use Poisson resummation 
to replace the sum by an integral and a sum over images </) ~ + Inr^p with p £ Z. The 
integral can again be done via contour deformation, replacing it by an infinite sum over 
residues as well. One then finds that the sum over the poles at a; = —im vanishes and we 
are left with the sum involving the quasinormal frequencies only. 



^ ^ ~ (,y.^i(2n+l+l±m)(t-iet)+im4> '^i'^ + n + l± m)T{l + U + I) 

7rr(/)r(/ + l)r+ 4^ 4^ ^ ^ ' r(l + n±m)r(l +n) 



where the ie factor is uniquely fixed by requesting convergence away from contact points 
and we suppressed the aforementioned sum over images. This expression can be evaluated 
without too much difficulty and adding the sum over images (remembering that the Poisson 
resummation yields an extra factor of r_|_), we finally get 

{TO{x)O{0)) = y . ,^^^^^^7). ^ • (4.4.43) 

^ ^ ^ " ^ [-cosh(t-iet) + cosh((?!. + 2m7rr+)]'+i ^ ^ 

This computation was done using the metric in (j4.4.4p where the mass of the BTZ 
entered through the periodicity of the angular coordinate (j4.4.5p . To restore standard 
conventions, we now perform the diffeomorphism t — > r+i and (/) ^ r^(j) followed by a 
Weyl transformation so that the boundary background metric is ds'^ = —dt^ + dcjp' with 
~ + 27r. Implementing these transformations in the two-point function we obtain 

(TO(x)O(O)) = y {2.Tr-HVi2^.) 

^1 ^_ ^^gj^(27rTt - iet) + cosh(2^r(,/. + 2m7r))]'+i ' ^ ^ ^ ^ 

where we reinstated the temperature T given in (j4.4.2|) . This correlator satisfies the KMS 
condition and is a sum over images in the (j) direction. It was obtained earlier via an analytic 
continuation of the Euclidean correlator in [38]. As discussed in more detail in [l2], it is 
related to the thermal AdS two-point function by a double analytic continuation. This 
can directly seen from (j4.4.43p . where the substitution t ^ icp and (j) ^ it yields precisely 
(j4.1.35p (up to ie insertions which then have to be inserted by hand). This is the real-time 
manifestion of the fact that Euclidean thermal AdSa and Euclidean BTZ, which are both 
filled tori, are related by an S transformation of the boundary torus. 
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Let us also write down the Wightman function, which can be obtained following the 
steps in section 14.1.21 

{0{x)0{x')) = (rcO[2](x)qi](x')) = ^^J] /'(ia;e--(*-*')+'^(^-*')x 

[a[2]-| a{uj, k, 1)(3{ljJ, k, I) + ap] a(— w, k, /)/?(— w, k, I)] . (4.4.45) 

We can again obtain a position-space expression by closing the contour and picking up the 
poles, which results in 

(O(,)O(0)) = V {2.Tr^HVi2^n) 

^^yx)^^u)j ^ ^_ ^^gj^^2^^^ _ .^^ ^ cosh(27rr((/) + 2m^))]'+i ' ^ ^ ^ ^ 

Finally, the retarded two-point function is of the form 

iAii{x, 0) = e{x){[0{x), 0(0)]) = {TO{x)O{0)) - {0{0)0{x)) . (4.4.47) 

From the above expressions, we find that it is analytic in the upper half of the complex 
w-plane and vanishes for t < 0. Actually, it has support only on the forward lightcone, 
which agrees with QFT expectations. Notice also that there is no need to insert ie's in the 
frequency-space expressions, since the poles in the complex frequency plane all have non- 
zero imaginary part. Such behavior however cannot arise in a CFT with a discrete energy 
spectrum, at least at finite N, where one expects that retarded correlators have poles on the 
real axis. Reconciling this behavior with expectations from the AdS/CFT correspondence 
is still an open issue; we refer to [111 [52] for discussions of this point. 

Let us finally remark that the retarded two-point function (j4.4.47p can also be shown 
to be related to purely ingoing boundary conditions at the horizon [13], leading eventually 
to the recipe of [12] . The derivation of this result from the current perspective is presented 
in more detail in [49] . 

4.5 Rotating black holes 

In the previous examples, we started with a CFT contour and obtained a corresponding 
bulk solution by the condition that it 'filled' this contour. In this section we will do the 
converse. We will start from a Lorentzian solution and look for Euclidean solutions that can 
be matched to it. This then leads to a specific CFT contour corresponding to the combined 
solution. 

Let us discuss the practical use of this procedure. As discussed in section [21 the parts 
of the solution associated with vertical segments of the contour are directly related to the 
initial and final state or density matrix of the field theory. The same information is also 
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encoded in the asymptotics of the Lorentzian solution, since from those one can compute 
the holographic 1-point functions and from them in principle one can reconstruct the dual 
state. The continuity of 1-point function across the matching surfaces guarantees that the 
information encoded in the Euclidean parts and the asymptotics of the Lorentzian solution 
is indeed the same. Typically, it is not very easy to extract the dual state starting from the 
vevs. The real-time methods discussed here present a new tool, namely given a Lorentzian 
solution one looks for Euclidean solutions that can be matched to it. One then uses this 
information to infer the holographic interpretation of the solution. 

In this subsection we illustrate how this is done using the rotating BTZ black hole 
[50\ [5T] . This discussion readily generalizes to higher dimensional rotating AdS-Kerr black 
holes [52l [53l [5l] . As one may expect, the contour turns out to be a thermal contour with a 
chemical potential for angular momentum. Furthermore, this example illustrates a number 
of additional issues as it provides a concrete example of the use of a complex metric. 

4.5.1 Lorentzian solution 

The metric for the three-dimensional rotating BTZ black hole \50\ [5T] is given by 



ds^ = -{r^ -rl- rl)dt^ + r^d^^ + 2r+r_dtd(l) + ^TT^ 2T ' (4-5-1) 



with (j) periodic, 

(j) (j) + 2tt . (4.5.2) 

The mass, angular momentum and temperature of the black hole are related to r+ and r_ 
via 

9 9 9 9 

It is convenient to use (t, (p, r) coordinates: 

i = r+t -I- r^4) , 
(j) = r^t + r+(^ , 

„2 „2 



(4.5.4) 



r2 



r — rz 



J.2 _ ^2 



Then the metric becomes 

df"^ 

ds"^ = -(f2 - l)dP + + f'^dct)^ , (4.5.5) 

— 1 

with the periodicity condition 

(t,(^) ~ (t-F27rr_,<^-h27rr+), (4.5.6) 
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with r_ and r_|_ real, and we consider < |r_| < but not the extremal case where 
|r_| = ?"+■ 

We consider an eternal rotating BTZ black hole with two radial boundaries. The rotat- 
ing BTZ black hole has a Penrose diagram that can be extended indefinitely to the future 
and the past, across the various horizons [51] • We will however cut off the spacetime along a 
spacelike hypersurface extending from one radial boundary to another, just as for the static 
BTZ example of the previous subsection. We thus explicitly avoid these extra regions and 
the singularities. 

4.5.2 Euclidean solution 

To find a boundary contour corresponding to this spacetime, we will first look for a Euclidean 
solution that is to be matched to the Lorentzian solution across some initial hypersurface. 
Usually, in passing to the Euclidean version of a rotating black hole we not only make 
the replacement t = —ir, but also analytically continue the angular momentum parameter 
(which is J or r_. in our case) to imaginary values. This way, the Euclidean metric one 
obtains is real. We will however show that the matching conditions are only satisfied for a 
complex Euclidean metric, given in coordinates (r, r, ip) by 

ds'^ = (r^ -rl- rl)dT^ + r^dif"^ - lir+r^drdip + _ 2 _ 2\ ' (4-5-7) 

with coordinate ranges we make precise below. We now discuss this metric in more detail. 

First of all, the Einstein equations are still satisfied for the complex metric, since they are 

satisfied for any real or complex r„. Second, a coordinate singularity arises at the horizons. 

Insisting on a nondegenerate metric, a (complex) coordinate transformation near the horizon 

shows that the necessary periodicity in Euclidean time that avoids such a singularity is 

27rr_i_ 27rir_ , 

{t,^)^{t+ ^ 2 '^^ + ^ ^' 4.5.8 

— r_ ''"+ ~ ^- 

which notably involves a translation in the imaginary 99 direction. To comply with this 
periodicity, we will take the Euclidean manifold Mq to be defined as follows. We first 
introduce Mq by extending the coordinates (r, r, y?) to complex values, with the periodicities 
as above. The metric (|4.5.7p should then be seen as a nondegenerate holomorphic (2, 0)- 
tensor on Mq. Within Mc, we take Mq to be the submanifold given by real r and r, but 
Im ip = T{r-/r^). Notice that Mq has three real dimensions. The metric restricts to Mq 
as a complex tensor and the volume element is a three-form which we can integrate along 
Mq. If we introduce 

IT 

if = 4) r, (4.5.9) 

r+ 
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then r, r and (p are real on Mq and therefore constitute an ordinary real coordinate system 
on Mq. In these coordinates the periodicity becomes 

(r,^)~(r + -^^,^) (4.5.10) 

and (p if! + 2-11 as well. However, the boundary metric in (r, ip) coordinates is no longer 
diagonal. Since this will complicate the analysis below, we continue to use the complex (p 
coordinate instead. 

4.5.3 Matching 

Let us now glue the 'Euclidean' and the Lorentzian manifolds together. We will first match 
the manifolds along a slice of constant t or r away from the horizon. Afterwards, we will 
deal with the subtleties introduced by the horizon. 

We begin with the first matching condition. On the Lorentzian side, we find that the 
induced metric on a slice of constant t is given by: 

hAsdx dx^ = r^d(f)'^ + — n^, n . (4.5.11) 

(r^ — rj^)[r^ — r±) 

On the Euclidean side, we find exactly the same metric on a slice of constant r, with 
the replacement — > 93, and therefore the first matching condition is satisfied. Let us 
now discuss the matching of the canonical momenta. On the Lorentzian side the extrinsic 
curvature (defined as usual with a real outward pointing unit normal) is 

^K,ABdx^dx^ = — =drd(l) (a<^'\9'\ 

and on the Euclidean side we obtain 

^ICABdx^dx^ = '^ir+r. ^^^^ ^ (4.5.13) 

/ (^2 _ j.'^'^^j.'i _ ^2 ^ 

A short computation then shows that the second matching condition is satisfied as well, 
including the factor of i. 

Notice that a complex metric is already needed in the first matching condition, i.e. the 
continuity equation for the induced metric Kab- Had we continued r_ to imaginary value 
on the Euclidean side so that the bulk metric is real, the induced metrics on the matching 
surface would not be the same (because the factor l/(r^ — r^) in (j4.5.11|) would become 
l/(r^ + r'^))- The fact that the metric is complex is therefore not directly related to the 
factor of i appearing in the matching conditions for the conjugate momenta. 
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Figure 13: The contour for the boundary CFT that corresponds to a rotating black hole 
does not only lie in the complex t plane, but also extends into the complex cj) plane. The 
circles should be identified. 

The matching conditions should also be satisfied at the horizons. Based on the example 
of the static BTZ black hole, one may consider using half a period of the Euclidean solution 
and matching the surface given byT = OtotL = and the surface given by r = 7rr+/(r^ — 
r^) to tji = 0. However, moving from r = to r = 7rr+/(r^ — r^) on Mq also involves 
an extra shift in the complex if direction. Therefore, a correct matching can be obtained 
by setting (j)^ = (p on the matching surface at L, and = (p + nir^/ir'^ — r^) at R. One 
may then verify that the matching condition are satisfied at the horizon by transforming 
to a coordinate system that is nonsingular at the horizon. Since d^p is a Killing vector, the 
matching conditions for gravity are insensitive to the extra twist in (p. However, for a scalar 
field the first matching condition becomes 



$o(r = 0, = 0, r) = $i,L(tL = 0, = 0, r) , 
^>o(t = 2 2 ' = -2 — ^) = '^iM^R = 0,(p = 0,r). 



(4.5.14) 



which is clearly sensitive to the twist in (p. 

We have shown that (half of) Mq with the metric (j4.5.7|) can be matched to the 
Lorentzian rotating BTZ black hole. We can therefore also glue two Euclidean and two 
Lorentzian spacetimes in the same manner as shown in Fig. [12] for the BTZ black hole. 
Analyzing then the boundary of this combination of spacetimes, we can finally read off the 
boundary contour corresponding to the rotating BTZ black hole: it is the contour of Fig. [T3j 
This is the same contour of Fig. [9l except that the vertical segments involve a shift in the 
imaginary (f) (or ip) direction of total magnitude 27rr_/(r^ — r^). Let us now interpret this 
result in field theory. 

First of all, notice that the boundary metric on the vertical segments is already complex, 
as can be verified by using real coordinates for the boundary of Mq via the coordinate 
transformation (j4.5.9p . This is in fact consistent with the anticipated result that this contour 
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corresponds to a CFT at finite temperature and with non-zero chemical potential for angular 
momentum. Namely, for such an ensemble the density matrix is 

p = exp{-f3{H + ^iP^)), (4.5.15) 

where H is the Hamiltonian and is a translation in (p. At the level of the path integral, 
such an ensemble corresponds to a contour that not only evolves in the imaginary time but 
also in the imaginary (j) direction. From the periodicity (j4.5.8p . we immediately read off: 

P=J^^ H = —. (4.5.16) 

— r_ r_|_ 

Of course, if one works purely in Euclidean time, one may also analytically continue ^u, r_ 
and J and then both the boundary and the bulk metric would be real. Our aim here was to 
develop a real-time formalism and this led to complex metrics both in the boundary theory 
and in the bulk spacetime. 

Let us finish this section with a brief comment on the use of complex (but non-degenerate) 
metrics in quantum gravity. First, it has been argued in the past (see for example [32]) 
that use of complex metrics after Wick rotation might be essential for a path integral over 
metrics. Second, saddle-point approximations often involve a deformation of the integration 
contour to a point in the complex plane, even if the integral originally is along the real axis. 
One particularly elementary example where this happens is a discretized vacuum-to- vacuum 
path integral for the harmonic oscillator, where the initial and final vacuum wave functions 
require such a contour deformation. Finally, complete reality of the bulk metric can no 
longer be maintained when one studies perturbations, as the ie insertions that follow from 
our prescription necessarily yield a complex graviton propagator. 

5 Summary of results and outlook 

We have presented a general prescription to holographically compute real-time correlation 
functions within the supergravity approximation. The main challenge in developing such 
a real-time prescription, relative to Euclidean methods, was to understand in detail how 
to deal with initial data. Our prescription is a direct 'holographic lift' of QFT real-time 
techniques to the gravitational setting, namely there is a gravitational counterpart of all 
QFT steps involved in such computations. In more detail, in QFT one typically chooses a 
contour in the complex time plane which usually consists of a sequence of horizontal (real) 
and vertical (imaginary) segments, the latter being related to the choice of density matrix or 
initial/final state. On the gravitational side, we construct solutions that directly correspond 
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to such QFT contours. Typically, real segments are associated with Lorentzian solutions and 
vertical segments with Euclidean solutions, with appropriate matching conditions imposed 
on the joining surface. The Euclidean parts encode the initial and final state in the field 
theory and this is refiected in the bulk, where they can be thought of as Hartle-Hawking 
wave functions. These wave functions also provide the necessary initial and final data for 
the perturbations around a given supergravity background. 

For the prescription to well defined, one must establish that one can remove all infinities 
through a process of (holographic) renormalization. Relative to the Euclidean discussion, 
new infinities can appear at timelike infinity. In our setup the analysis boils down to 
analyzing possible new contributions from the joining surfaces. We show that no new 
counterterms are needed and the holographic 1-point functions are continuous across the 
matching surface. The continuity of the 1-point functions is an important consistency 
condition of the entire setup: as mentioned above, the Euclidean parts of the solution are 
directly related to the initial/final state but as is also well known the 1-point functions 
encode the same information, too. 

As a sidenote, the holographic nature of the prescription also nicely shows up in the 
following issue that we encountered when demonstrating the renormalization. Starting from 
a boundary state defined at a boundary Cauchy surface, say the surface t = to, one can 
extend this surface to the bulk, t = f{r, x) with /(r, x) ^ to as r ^ oo, but clearly there is 
a certain amount of freedom of how this is done, parametrized by the subleading behavior 
of /(r, x). These extensions are not part of the boundary theory, so the renormalized theory 
should be independent of them. We explicitly find that possible dependence on /(r, x) drops 
out indeed. 

Having set up the prescription, we then moved on to demonstrate how to apply it in 
a variety of examples that each illustrate different points. The first example involved the 
holographic computation of a vacuum Wightman function. Although its functional form 
was already known, we were able to compute it completely holographically, without analytic 
continuation or insertion of ie by hand; instead, our computation provided for all the right 
signs and ie insertions. We then computed a two-point function in real-time thermal AdS. In 
this computation, we used the same Lorentzian background but different initial data, which 
highlights the importance of properly defining the initial and final boundary conditions. 

The prescription can be used to compute higher-point functions as well and we explic- 
itly demonstrated how to do such computations in an AdS background. This discussion 
straightforwardly extends to any other Asymptotically (locally) AdS bulk spacetime. It is 
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worth mentioning that our prescription resulted in a bulk-bulk propagator which is already 
of quantum-mechanical nature {i.e. Feynman rather than retarded). This shows that the 
bulk fields are path-integral quantized and the Euclidean caps provide the proper initial 
and final states. The prescription thus naturally incorporates QFT in curved space and it 
is not necessary to quantize perturbations by hand again. 

A real-time thermal contour can also be 'filled' with an eternal black hole spacetime. 
Despite the presence of singularities and horizons, we demonstrated how the initial and 
final conditions could again be unambiguously specified via Euclidean caps. This procedure 
extends to rotating black holes, where the analytic continuation is more subtle. In our 
case, the reality conditions of the bulk fields and factors of i that arise in passing from real 
to imaginary time agree with QFT arguments, where the situation is well- understood. In 
particular, this procedure led to a complex bulk (and boundary) metric in the case of a 
rotating black hole. 

The correlators we computed in the various examples were largely known from earlier 
work, where they were obtained using special properties of the backgrounds and analytic 
continuation. The emphasis here was on the coherent derivation of these results using 
the new real-time prescription: statistical factors and appropriate ie insertions in 2-point 
functions all follow uniquely from solving the matching conditions. 

The true power of the new method however should be in the applications that lie ahead 
of us. Current and future applications of holography to RHIC and LHC physics or to 
condensed matter systems require holographic modeling of non-equilibrium phenomena and 
for such applications previous methods are just not applicable. On a more fundamental 
level, the new prescription may help us addressing global issues and questions regarding the 
holographic encoding of the bulk causal structure, including bulk horizons, and the parts of 
spacetime beyond the horizon. It would also be particularly interesting to extend the black 
hole analysis of this paper to a collapsing shell of matter aiming at a holographic description 
of the process of black hole formation. Work about some of these issues is under way. 
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A Real-time quantum field theory 



In this appendix, we discuss some aspects of real-time quantum field theory relevant for our 
discussion. The material presented here is not new and it is included to make this paper 
self-contained. 

A.l Vacuum wave function insertions 

In this section, we will analyze how the vacuum wave function insertions in the path integral 
lead to ie insertions. In the main text, we mentioned how the wave functions can be obtained 
as path integrals along vertical segments in the complex time plane, leading ultimately to 
a contour as in Fig. [1^. Let us begin by an explicit computation of these wave functions in 
a relatively simple case. 

A.1.1 Computation of the wave functions 

We will take a real free massive boson on flat Minkowski space M^'"^"^. As explained in the 
main text, the initial wave function ((/>_, —T\Q) is computed via the projection: 

lim e'^^""=((/._, -Tle-^^l^-) = -r|f])(!^|^) . (A.1.1) 

For simplicity, we shift the time coordinate such that —T — > and we will take |^) = {(pp, i[3) 
for some spatial field configuration (j)p{x). 

Since we take the field to be free, the path integral is Gaussian and can be computed 
exactly. Let (j){t,x) be the solution to the equation of motion satisfying <j){if3,x) = (t)p{x) 
and (/>(0, x) = (/>-(x). We then obtain 

(0_,O|O) = lim A^e-^^[^1 , (A.1.2) 

with a normalization that does not depend on (j)^ and Se the Euclidean on-shell action 
for the boson. Introducing a Euclidean time coordinate r = it, this action is given by: 

Se[(I>\ = \f<^^j '^'^'^^ id,.<pd^'(t> + m^cl?) . (A.1.3) 
On-shell, it reduces to a surface integral, 

0|J]) = lim A^exp (- \ f (f-^x [0(r, x)a.0(r, x)f.) . (A.1.4) 

Finding (j) is not hard and in the limit /3 — > oo we find that all dependence on cpp can be 
absorbed in a shift of and we recover the usual Gaussian wave function (5^, written in 
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Fourier space as 

((A_,0|^^)=AA'exp(-i I ^^1^0_(fcK</>„(-A;)), (A.1.5) 

with u>k = Vfc^ + ?7T.2 . The conjugate final wave function {^}\(f)^,T) can be computed using 
the same procedure, leading to exactly the same result. 

If interactions are switched on, the wave functions receive corrections. However, as long 
as these interactions can be switched off adiabatically for large times, the corrections can 
also be ignored in the limit oo. The analogous case in thermal field theory, which 

we discuss below, is briefly discussed in \A7\ section 2.4.1]. For massless field theories there 
are subtleties, but these considerations are not directly relevant for us and they will not be 
discussed here. A computation of the ground state wave function for electromagnetism and 
linearized gravity can be found in [56j. 

A. 1.2 Effect of the wave function insertions 

Let us now show how the wave function insertions determine ie-insertions in the propagator. 
To this end, we introduce a source J and compute 

Z[J] = {n\e-'^'^'''\n) . (A.1.6) 

We suppose that the source vanishes smoothly at the endpoints t = itT of the Lorentzian 
segment. Again via the usual slicing arguments, the path- integral representation one obtains 
is 

ZIJ] = j [V^] exp - i I dtd'^-^x J4>-1[ -^-^<t^-{k)u:k'i>-{-k) 



1 f d'^-^k 



2 



P"'' (A.1.7) 



where (j)±{k) is the Fourier transform of 0(ibT, x) with respect to the spatial coordinates. 
Notice that the boundary values for the path integral J['D(j)] are not fixed. 

To compute the path integral, we shift the integrand (f) = X + V'? where x satisfies 
□x — m?x = J cind tp is the new integration variable. Notice that x is not uniquely defined 
unless we specify some boundary conditions. To find these, notice that the aim of this shift 
is to get all the factors involving J and x to come out in front of the path integral, resulting 
in 

Z[J] =AAexp( - ^ I d'^xxj) , (A.1.8) 

from which we would directly obtain the propagator as is shown below. However, an 
analysis of the boundary terms shows that such a factorization only occurs if one imposes 
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additionally the two extra constraints: 



(A.1.9) 



(2vr) 

which should hold for all values of 'i/'i- These conditions provide the boundary conditions 
for X- Since the source vanishes at the endpoints, x is homogeneous for t = itT, and 
therefore has a Fourier expansion involving only modes of the form ^T^'^kt+^kx ^ The boundary 
conditions that one derives from these constraints are then simply that xi~T,x) should 
contain only negative frequencies (i.e. modes of the form q-'^^^+^^x ^j^j^ u; < 0) and x{T, x) 
should contain only positive frequencies. But this uniquely fixes x to be of the form 

X = j dt'(f~^x'I\F{t-t',x-x')J{x'), (A.1.10) 

with /S.F the Feynman propagator, 

, , , f dtd'^-^x e-i^t+ikx 

AF(t,x) = / — — -j ^ ^ . (A.1.11) 

As one may verify by contour deformation, one indeed obtains only positive/negative fre- 
quencies to the future/past of the source. 

We can now take the limit T ^ oo. Assuming that the source and any perturbatively 
added interactions vanish slowly at late times, the propagator and the wave functions are 
unmodified and all that is left are the ie-insertions which enter in the perturbative expansion, 
which is precisely what we wanted to show. 

Different (equivalent) arguments that translate wave functions to ie insertions can be 
found in the textbooks [55] and [57] . In particular in |57j , the contour of Fig. [1^ is deformed 
to a straight line that runs almost parallel to the real time axis, from — T(l — ie) to T(l — ie), 
with T ^ oo. The projection property is left unchanged and this way one still obtains 
vacuum-to- vacuum amplitudes. The contour should always go downward or horizontal in 
the complex time plane so that the operator exp(—iH At) remains finite. 

Finally, notice that the saddle-point x is actually a complex solution, although we started 
with a real scalar field and a real source J{x). This can be viewed as a contour deformation in 
field space before taking the saddle-point approximation. Such a deformation is very explicit 
when one discretizes the path integral. Nevertheless, the usual hermiticity constraints of 
n-point functions are still satisfied. The fact that a saddle-point approximation may involve 
complex fields holds for gravity as well. In the context of holography, it is the hermiticity 
of the boundary stress energy tensor and its correlators that restricts the allowed complex 
metrics. 
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A. 2 Contour time ordering 

The prescription above showed that time-ordered vacuum-to-vacuum correlation functions 
can be obtained via a path integral along a specific contour in the complex time plane. If 
the initial state is not an energy eigenstate then a corresponding path-integral formulation 
involves a so-called in-in contour which runs from, say t = to t = T and back to t = 0. 
At the endpoints of the contour, one may impose initial and final conditions on the fields 
corresponding to the state or ensemble under consideration. 
As an example, suppose one wants to compute 

{^\0{t')\^) = (^|e*^*'Oe-*^*'|^) . (A.2.1) 

We see that we first have to evolve the state |^) for a time t' before we insert the operator, 
but afterwards we also have to evolve back in time before we insert the final wave function. 
This is the in-in or 'closed time path' formalism of |27 t l28 t[29ll30j . Extending the contour to 
go beyond the point t', say to some point t", and then back again amounts to an insertion 
of the identity in the form exp[iH{t" — t')]exp[—iH{t" — t')]. Such an extension of the 
contour will not change the overall amplitude, which is something that is reflected in the 
dual gravity theory as well. 

Another example where the in-in formalism is useful is real-time quantum field theory 
at finite temperature. In that case, the ensemble is described by a thermal density matrix, 

p = exp{-(3H) , (A.2.2) 

with (3 = 1/T and H the Hamiltonian. Expectation values in such an ensemble are traces: 

{0)f3 = Tr{e~^^0). (A.2.3) 

One-point functions can be computed using the path-integral formalism by taking a contour 
that runs straight down along the imaginary time axis from to with (anti-)periodic 
boundary conditions for the fields. However, a convenient way to obtain dynamical informa- 
tion, i.e. real-time correlation functions, is again via the in-in formalism. In that case, one 
still rewrites the density matrix as a Euclidean path integral, but the operator insertions 
at different times force one to also path integrate along an in-in contour running along the 
real time axis. The resulting contour is drawn in Fig. [TJ; on page[9l 

A.2.1 Two-point functions 

For an in-in contour, one may insert operators along both the forward- and the backward- 
going segments of the contour. Via the usual slicing arguments, the correlation functions 
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from a path integral along this contour are contour-time- ordered. That is, if we pick a real 
'time' parameter tc that increases monotonically along the contour, then the correlation 
functions obtained from a path integral along this contour are ordered from small to large 
tc- For example, let us take an in-in contour as in Fig. [T]3. The Lorentzian part runs up 
and then down the real t axis, from to T and then back to 0. We denote the first segment 
by Ci and the second by C2. One may then choose tc = ^ on Ci and tc = 2T — t on C2. 

Let us now introduce a source J along the contour to compute the two-point functions. 
For an in-in contour, the source-operator coupling —i dtJO can be split in two parts and 
the partition function is defined as: 

^[•/[ih^p]] = (exp(-i^ dtid'^-^xJ[^0[^]+i dt2d'^-^xJ[2]0[2])) , (A.2.4) 

where a subscript in square brackets denotes the segment on which the field lives. The 
expectation values of course depend on the ensemble or state that is specified at t = 0, but 
we will not write this explicitly. 

Via functional differentiation one obtains four possible two-point functions, 

with Tc denoting contour-time-ordering. Along the first segment contour-time-ordering 
coincides with normal time-ordering, 

(TeO[i](x)0[i](x')) = (TO(x)O(x')) . (A.2.6) 

Along the second, backward-running segment, contour-time-ordering coincides with anti- 
time ordering, denoted by T, 

(T,0[2](x)0[2](x')) = {fO{x)0{x')) . (A.2.7) 

If one puts one argument on the forward contour and the other on the backward contour, 
the latter one will always be later in contour time than the former and we get the Wightman 
functions: 

(T,0[i](x)0[2](x')) = (0[2](x')0[i](x)) = {0{x')0{x)) , ^^^^^ 
(rcO[2](x)O[i](x0) = (0[2](x)0[i](x')) = {0{x)0{x')) . 

Notice that the in-in path is also suitable to obtain vacuum-to-vacuum Wightman func- 
tions from a path integral. In the main text, we perform the corresponding holographic 
computation. 
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A. 2. 2 Linear response 

Finally, we discuss the important role of the retarded two-point function, which describes 
the reaction of the system to an external perturbation. The perturbation can be described 
as a deformation of the theory such that the action S changes to S — J JO. In the in- 
in formalism, the deformation should be present on both contours, so J[i] = Jp] = J- 
Expanding then in J, we obtain the first-order response to the one-point function on Ci: 

5{Oix'))j = J^dxJ{x)j^^{0[,]{x'))j 

= -i r dhJ{x){%0[^^{x)0[^ix'))+i r dt2J{x){%0[2]ix)0[^{x')) 

Jo Jo 

r-T rO 

= -i dtJ{x){TO{x)0{x')) +i dt{0{x)0{x')) 

Jo JT 



(A.2.9) 



T 



dt J{x)Afi{x' , x) 







where we suppressed the spatial integrations and Aji{x,x') is the retarded two-point func- 
tion, 

Ar{x',x) = -i9{x' -x){[0{x'),0{x)]) , (A.2.10) 
which vanishes outside the future lightcone. The response is thus causal, as expected. 

A. 3 ie-insertions 

In this section, we briefly review the analytic properties of two-point functions and corre- 
sponding ie-insertions. 

We start with the Wightman function 

(^(x)^(O)), (A.3.1) 

which is analytic in the lower half of the complex t plane [58] . The Wightman function can 
be obtained by the replacement —it = t — ie in the Euclidean correlator, because then the 
poles along the real t axis are shifted into the upper half of the complex t plane. Its Fourier 
transform, 

'"dtd'^-ixe*"*-*^'^-(V'(^)V'(0)) , (A.3.2) 



vanishes for negative frequencies, since we can close the contour for the t-integral via the 
lower half plane. Positivity of the spectral density also implies that the Fourier transform 
is a real and positive distribution for positive frequencies [58]. The Fourier transform thus 
maps a function (or distribution) that is analytic in a upper or lower half plane to a function 
that vanishes on the right or the left real axis. 
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Next, the time-ordered two-point function is defined as 

(rV(x)V(O)) = 0(t)(V(x)V(O)) + 0(-t)(V(O)V(x)) , (A.3.3) 

wfiich can be obtained from the Euchdean correlator by the replacement —ir = t — iet. Its 
poles are shifted into the upper half of the complex t plane for Re t > and in the lower half 
plane for Re t < 0. To obtain the Fourier transform, we close the contour in the appropriate 
half plane in the complex time plane. Picking up the poles, we find a sum over positive 
frequencies for t > and one over negative frequencies for t < 0. This implies that we 
need the usual Feynman contour around the poles to define the inverse Fourier transform. 
One may replace to ^ uj + ieto in the Fourier-space expression to explicitly indicate such a 
contour. Obviously 

- {lo + ietof = -io'^ - ie (A.3.4) 

and for example the propagator (jA.l.llh indeed has the required analyticity properties. 

The retarded two-point function was defined in (|A.2.10p . Causality of the field theory 
determines that it vanishes completely outside the future lightcone. We may write it as an 
inverse Fourier transform: 

Ar{x, 0) = ^ J dc^d^-ifc e-*"*+*^^Afi(a;, k) . (A.3.5) 

Notice that A/{(x,0) vanishes for t < 0. Since we can then close the (j integral in (|A.3.5P 
in the upper half plane, we find that Af{{uj,k) must be analytic in the upper half of the 
complex frequency plane. Finally, the advanced two-point function is the reversed retarded 
two-point function: 

Aa{x,x') = Ar{x',x). (A.3.6) 

It therefore vanishes outside of the past lightcone and is analytic in the lower half of the 
complex frequency plane. 
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